AFAMRL-TR-85-015 


AD-A272  919 

lllllli 


SPACING  POINTS  IN  A  THREE  DIMENSIONAL  CONVEX  REGION 
FOR  MAXIMUM  SEPARATION:  A  COLOR-SPACE  APPLICATION 


ROSSE.  ROLEY,  CART 


AIR  FORCE  AEROSPACE  MEDICAL  RESEARCH  LABORATORY 
Wright-Patterson  Air  Force  Base,  Ohio  45433 


APRIL  1985 


Approved  for  public  release;  distribution  unlimited. 


AIR  FORCE  AEROSPA CE  MEDICAL  RESEARCH  LABORATORY 
AEROSPACE  MEDICAL  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHTPATTERSON  AIR  FORCE  BASE.  OHIO  45433 


93-28483 


NOTICES 


When  US  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  a  definitely  related  Government  procurement  opera¬ 
tion,  the  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever,  and  the  fact  that  the  Government  may  have  formulated,  furnish 
ed,  or  in  any  way  supplied  the  said  drawings,  specifications,  or  other 
data,  is  not  to  be  regarded  by  implication  or  otherwise,  as  in  any  manner 
licensing  the  holder  or  any  other  person  or  corporation,  or  conveying  any 
rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention 
that  may  in  any  way  be  related  thereto. 


Please  do  not  request  copies  of  this  report  from  Air  Force  Aerospace  Med¬ 
ical  Research  Laboratory.  Additional  copies  may  be  purchased  from: 

National  Technical  Information  Service 
5285  Port  Royal  Road 
Springfield,  Virginia  22161 

Federal  Government  agencies  and  their  contractors  registered  with  Defense 
Technical  Information  Center  should  direct  requests  for  copies  of  this 
report  to: 


Defense  Technical  Information  Center 
Cameron  Station 
Alexandria,  Virginia  22314 


TECHNICAL  REVIEW  AND  APPROVAL 

AFAMRL-TR-85-015 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS 
it  will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 

FOR  THE  COWANDER 


CHARCB  BATE! . 

Director,  Human  Engineering  Division 

Air  Force  Aerospace  Medical  Research  Laboratory 


SfCuHiTv  CvASSii'iCA'^ON  Of  TmiS  gage 


1*.  «ii»o«T  secuniTv  classiP'Cation 

Unclas’^ified 


2a  seCoAiTv  CUASSl^ 'CA  r  lON  AvjTmOBiTV 


I 


REPORT  DOCUMENTATION  PAGE 


lb.  AESrniCTlvC  MAAKlNCS 

None 


2b.  OeCLASSlP'CATiON  OOWVNOAAOiNC  SCMEDULE 


4  PERsonMiNG  oRGaniZat.on  REPORT  NoMBERiSI 

AFAMRL-TR-85-015 


6«.  name  op  PERPORMiNO  organization  So.  oppiCE  SvmbOl 

Air  Force  Aerospace  Medical  ifapphc^bit, 

Research  Laboratory 


6c.  AOOAESS  'City  S(a<tf  o/id  /.IP  CMi4t 

Wri ght-Patterson  Air  Force  Base 
Ohio  45433 


NAMC  FuNOiNG^SPONSOniNC 
organization 


86  OPPiC£  Sv  MSOl 
tif  Afipltcabist 


8c.  aOOR£S$  (Cfy  :3tat€  and  /IP  Ca<Ut 

Wright-Patterson  Air  Force  Base 
Ohio  45433 


3  O-STRieoTION/AVAlLABlLiTV  OP  REPORT 

Approved  for  public  release;  distribution  is 
uni imited. 


6.  monitoring  organization  report  NoMSERtS) 

N/A 


7a.  NAME  OP  MONITORING  ORGANIZATION 

N/A 


7b.  address  (City.  Stats  and  ZIP  CodSf 

N/A 


9.  PROCUREMENT  INSTRUMENT  lOENTiPiCATION  NUMBER 

F33615-82-C-0511 


to  SOURCE  OP  PUNOING  NOS 


program 

FROjECT  1 

task 

MORK  UNIT 

Element  no. 

NO. 

NO. 

NO. 

62202F 

7184 

11  1 

! 

49 

12.  personal  auTmQRiSI 


laniawiii 


13d  Tvpe  Of  REPORT 


I  ■uFnirniinBliKCEnkl 


16  SUPPLEMENTAR  Y  notation 

None 


17  COSATI  COOES 


t.  USAF 


130,  time  COVERED 

PROM  9/84 


12/84 


14  DATE  OP  REPORT /Tr  .  .Wo..  0«yi  IS.  PAQE  COUNT 

1985  April  76 


18.  SUBJECT  TERMS  JContmus  on  if  ntesMSQry  and  idsnttfy  l>y  Mocd  numba^t 

Algorithm  Color  Selection  Display  Maximin  Spacing 

Color  CRT  Distance  Maximum 

f.n To r  Coding  DISC r iminat i 0 n  Heuristic  Optimal _ 


19.  abstract  Con/mutf  on  f  ntct$»ary  and  <d«nn^y  by  btoeb  numbari 

This  report  introduces  a  new  method  for  solving  the  problem  of  optimally  spacing  points  in 
a  three-dimensional  region  so  that  their  distances  from  each  other  are  as  great  as  pos- 
sible.  One  application  of  the  problem  deals  with  color  selection  for  aircraft  displays 
where  the  colors  are  plotted  as  points  in  a  three-dimensional  color  space  and  the  distance 
between  two  points  is  directly  related  to  the  distinguishability  of  the  two  colors.  The 
method  itself  is  a  heuristic  algorithm  very  similar  to  one  designed  by  Carter  and  Carter. 
The  newer  algorithm  apparently  yields  similar  solutions  with  fewer  runs,  but  because  it  is 
more  thorough,  it  is  slower.  The  program  was  tested  on  problems  as  large  as  23  points 
whose  feasible  region  had  seven  faces.  The  major  disadvantage  of  this  new  method  is  that, 
like  Carter  and  Carter's,  its  solutions  are  not  guaranteed  to  be  optimal.  As  a  result,  the 
user  must  still  perform  several  replications  at  various  randomly  selected  starting  loca¬ 
tions  in  order  to  increase  the  chances  of  achieving  an  optimal  solution. 


2a  oistribution/availability  of  abstract 

UNCLASSIFIEO/UNLIMITEO  □  SAME  AS  RFT  ^OTlC  USERS  G 

21.  abstract  security  classification 

Unclassi fied 

23a  name  of  RESFONSiBlE  INGIVIOUA.. 

23b.  TBLEFMONS  number 
i/ncluW*  Art*  Coda) 

32c.  office  symbol 

David  L.  Post 

(5131  255-7597 

AFAMRL/HEA 

DO  FORM  U73,  83  APR 


EDITION  OR  1  JAN  73  IS  OBSOLITI. 


SECURITY  CLASSIFICATION  OF  THIS  FACS 


SieuAITV  CLAa«l»iCATION  Of  THIS  >AO« 


Block  11  (continued) 

SPACING  POINTS  IN  A  THREE-DIMENSIONAL  CONVEX  REGION 
FOR  MAXIMUM  SEPARATION:  A  COLOR-SPACE  APPLICATION 


PREFACE 


This  report  describes  a  project  that  was  performed  in  partial  ful  f i  I  Imetit  of 
requirements  for  the  degree  of  Master's  of  Philosophy  in  Uperations 
at  the  Air  Force  Institute  of  Technology,  Wright-Patterson  Air  Force  Ba-.e, 
Ohio  45433.  The  project  was  sponsored  by  the  Visual  Display  Systems  Branch, 
Human  Engineering  Division,  Air  Force  Aerospace  Medical  Research  Laboratory 
(AFAMRL/HEA) ,  under  Project  7184,  Man-Machine  Integration  Technology,  work 
Unit  7184-11-49,  Color  Display  Laboratory. 

The  author  wishes  to  thank  his  thesis  advisor,  Lt  Col  Ivey  Cook,  and 
Dr.  David  L.  Post  (AFAMRL/HEA)  for  their  assistance  with  the  project. 


1 


TABLE  OF  CONTENTS 


Section  Page 

1  INTRODUCTION  5 

2  BACKGROUND  6 

3  LITERATURE  REVIEW  OF  CARTER  AND  CARTER'S  METHOD  11 

4  MATHEMATICAL  DEVELOPMENT  14 

PROBLEM  FORMULATION  14 

SOLUTIONS  TO  SIMPLE  PROBLEMS  22 

SOLUTIONS  BY  NONLINEAR  PROGRAMMING  26 

5  DESCRIPTION  OF  HEURISTIC  28 

THE  ALGORITHM  28 

DIFFERENCES  WITH  CARTER  AND  CARTER'S  METHOD  30 

6  RESULTS  33 

OPTIMALITY  33 

SIZE  LIMITATIONS  33 

PERFORMANCE  IN  COMPARISON  TO  CARTER 
AND  CARTER'S  METHOD  34 

SENSITIVITY  TO  PARAMETERS  38 

ADVANTAGES  40 

DISADVANTAGES  41 

7  SUMMARY  42 

CONCLUSIONS  42 

SUGGESTIONS  FOR  FURTHER  RESEARCH  43 

Appendix  A  COLOR  TRANSFORMATIONS  45 

Appendix  B  FORTRAN  CODE  AND  DOCUMENTATION  47 

Appendix  C  EXAMPLE  OF  ROLEY  HEURISTIC  SPACING  FOUR 

POINTS  IN  A  SQUARE  67 

BIBLIOGRAPHY  75 


2 


LIST  OF  ILLUSTRATIONS 


Page 

8 


Figure 

1  Example  Feasible  Region  for  the  Color-Spacing  Problem 

2  1976  CIE  UCS  Diagram  9 

3  Alternative  Locations  for  an  Endpoint  11 

4  Solution  for  Spacing  Two  Points  in  a  Square  23 

5  Solution  for  Spacing  Three  Points  in  a  Square  23 

6  Solution  for  Spacing  Four  Points  in  a  Square  24 

7  Solution  for  Spacing  Five  Points  in  a  Square  24 

8  Solution  for  Spacing  Six  Points  in  a  Square  25 

9  Approximate  Feasible  Region  for  the  Color-Spacing 

Problem  in  the  CIE  L*u*v*  Color  System  29 

10  Flowchart  of  Roley  Heuristic  31 

11  Example  of  a  Locally  Optimal  Solution  40 


3 


LIST  OF  TABLES 


Number 

1  COMPARISON  OF  MINIMUM  DISTANCES 

2  COMPARISON  OF  ALL  DISTANCES  FOR  THE  SIX-COLOR  PROBLEM 


4 


Section  1 
INTRODUCTION 


The  problem  addressed  herein  is  to  space  points  in  a  three-dimensional 
region  so  that  they  are  as  far  apart  as  possible.  A  simple  example  would  be 
to  maximize  the  distance  between  two  points  in  a  cube.  For  this  case,  the 
solution  should  be  fairly  obvious  (place  the  points  in  opposite  corners) 
but,  in  general,  solutions  to  this  spacing  problem  are  not  obtained  as 
readily  as  one  might  expect.  Consider,  for  example,  the  difficulties  that 
would  arise  if  three  points  must  be  spaced  in  the  cube.  If  the  region  were 
shaped  less  regularly  (e.g.,  like  a  septehedron),  the  problem  would  become 
even  more  complex. 

The  spacing  problem  that  is  explored  in  this  report  is  new  to  operations 
research.  It  was  introduced  by  Carter  and  Carter  (1982)  in  the  context  of 
color  displays,  but  has  received  no  other  attention  in  published  litera¬ 
ture.  The  relevance  to  color  displays  stems  from  the  fact  that  colors  can 
be  represented  as  points  in  a  three-dimensional  region.  If  the  spatial  dis¬ 
tribution  of  color  within  this  region  is  perceptually  uniform,  the  distance 
between  any  pair  of  colors  can  be  taken  to  represent  their  perceptual  dif¬ 
ference.  Thus,  if  one  wishes  to  select  several  maximally  discriminable 
colors  for  use  on  a  display,  one  can  approach  the  problem  by  attempting  to 
space  the  points  representing  those  colors  in  the  perceptually  uniform 
region  so  that  they  are  maximally  distant  from  one  another. 

The  present  report  also  discusses  the  spacing  problem  in  the  context  of 
color  selection.  It  is  hoped  that  this  will  aid  in  the  design  of  color  dis¬ 
plays  and  benefit;  for  example,  pilots  who  rely  on  color  to  quickly  distin¬ 
guish  between  hostile  versus  friendly  forces,  safe  versus  dangerous  flight 
conditions,  etc.  It  is  also  hoped  that  this  initial  exploration  will  con¬ 
tribute  substantially  to  the  understanding  of  the  general  spacing  problem 
and  thereby  lead  to  further  work  on  this  new  problem  in  operations  research. 


5 


Section  2 
BACKGROUND 


Why  would  anyone  want  to  space  points  in  a  three-dimensional  region?  There 
are  many  possible  reasons.  Perhaps  you  want  to  locate  speakers  for  a 
quadraphonic  stereo  system  in  your  living  room.  You  might  want  them  as  far 
apart  from  each  other  as  possible  to  maximize  the  size  of  the  soundstage. 

Or  perhaps  you  own  a  warehouse  with  a  limited  number  of  security  cameras. 

You  would  probably  want  the  cameras  spaced  so  that  no  two  were  filming  the 
same  areas  of  the  warehouse.  How  would  you  go  about  solving  this  problem 
mathematically?  You  could  try  to  maximize  the  sum  of  the  distances  between 
the  objects,  but  that  might  result  in  a  solution  with  an  unacceptably  small 
distance.  A  more  desirable  solution  would  probably  result  from  maximizing 
the  distance  between  the  two  closest  objects.  This  "maximin"  formulation  is 
accepted  as  the  objective  for  the  spacing  problem. 

The  spacing  problem  is  very  similar  to  a  problem  known  in  the  world  of  oper¬ 
ations  research  as  the  obnoxious-facility  location  problem.  An  example  of 
an  obnoxious  facility  would  be  a  nuclear  waste  disposal  site  which  is  most 
desirable  when  it  is  located  far  away  from  cities.  Church  and  Garfinkel 
(1978)  offer  a  solution  to  the  obnoxious-facility  location  problem  where  the 
facilities  can  be  placed  in  a  discrete  number  of  locations  on  a  network,  but 
the  differences  between  that  problem  and  the  spacing  problem  preclude  suc¬ 
cessfully  adapting  Church  and  Garfinkel 's  method  to  the  spacing  problem. 
First  of  all,  the  idea  is  to  locate  obnoxious  facilities  as  far  away  as 
possible  from  other  fixed  points,  whereas  the  spacing  problem  aims  to  locate 
the  objects  as  far  away  from  each  other  as  possible.  Second,  Church  and 
Garfinkel 's  facilities  can  be  placed  in  only  a  finite  number  of  locations  on 
a  network,  while  the  objects  in  the  spacing  problem  can  take  on  an  infinit<» 
number  of  locations.  Finally,  the  spacing  problem's  objective  is  to  maxi¬ 
mize  the  minimum  distance  between  points  (maximin),  as  opposed  to  Church  and 
Garfinkel 's  objective  of  maximizing  the  median  distance  (maxian). 

Another  problem  similar  to  the  spacing  problem  is  the  location  problem.  It 
seeks  to  locate  objects  (like  warehouses)  as  close  to  other  fixed  facilities 
(such  as  retail  stores)  as  possible.  Once  again,  much  work  has  been  done  in 
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the  area  by  many  people,  most  notably  Charalambous  (1981),  Calamai  and 
Charalambous  (1980),  Cooper  (1963),  Juel  (1982),  Love  (1976),  Love  and  duel 
(1982),  and  Love  and  Morris  '1975),  but  any  applicability  to  the  spacing 
problem  is  negligible.  As  uefore,  there  is  the  problem  of  locating  facili¬ 
ties  in  relation  to  fixed  objects  instead  of  to  each  other.  The  location 
problem  is  also  concerned  with  minimizing  the  sum  of  the  distances  (minisum) 
or  minimizi'^g  the  maximum  distance  (minimax),  not  with  maximizing  the 
minimum  uistance. 

There  are  many  practical  applications  of  the  spacing  problem.  In  addition 
to  those  mentioned  before,  the  spacing  problem  could  be  applied  toward 
locating  MX  missile  silos,  spacing  mines  in  a  mine  field,  or  perhaps  placing 
communications  satellites  in  space  for  maximum  coverage  of  the  earth.  A 
particularly  fascinating  application  of  the  spacing  problem  is  in  the  area 
of  color  spacing.  One  of  the  tasks  of  the  Air  Force  Aerospace  Medical 
Research  Laboratory  (AFAMRL)  is  to  choose  colors  for  aircraft  displays,  air 
traffic  control  displays,  aeronautical  maps,  etc.,  so  that  all  the  colors 
are  as  distinguishable  from  each  other  as  possible.  When  one  is  dealing 
with  a  color  cathode-ray  tube  (CRT)  display,  the  CRT’s  red,  green,  and  blue 
guns  are  used  additively  to  produce  various  colors.  Therefore,  any  color 
can  be  uniquely  defined  by  three  parameters.  They  are  the  luminances  of  the 
CRT's  red,  green,  and  blue  guns.  So,  each  color  produced  by  the  CRT  can  be 
plotted  as  a  point  in  three  dimensions  where  the  axes  represent  those  three 
luminance  values.  A  depiction  of  the  resulting  color  space  is  shown  in 
Figure  1.  The  boundaries  of  the  region  correspond  to  the  technological 
limits  of  the  CRT.  For  instance,  in  Figure  1,  the  maximum  luminances  are 
50,  160,  ano  20  candelas  per  meter  squared  for  the  red,  green,  and  blue 
guns,  respectively,  with  zero  being  the  minimum.  In  addition,  there  is  a 
small  region  near  the  origin  that  is  not  feasible  because  the  colors  in  this 
area  do  not  have  sufficient  luminance  to  be  readily  seen. 

Unfortunately,  this  red,  green,  blue  CRT-luminance  space  is  not  perceptually 
uniform.  That  is  to  say  that  the  perceived  difference  between  the  red  and 
the  purple  shown  in  Figure  1  is  not  the  same  as  the  difference  between 
yellow  and  white  even  though  their  distances  are  the  same.  The  Commission 
Internationale  de  I'Eclairage  (CIE)  recommends  the  1976  CIE  L*u*v*  system 


7 


Red 


Figure  1.  Example  Feasible  Region  for  the  Color  Problem 

as  a  more  perceptually  uniform  color  space  (Robertson,  1977).  The  L*  axis 
is  a  function  of  the  luminance  of  a  color,  while  the  u*v*  plane  identifies 
the  color's  position  in  a  transformation  of  the  1976  CIE  uniform 
chromaticity-scale  (UCS)  diagram  (see  Figure  2).  The  distinguishability  of 
two  colors  in  the  CIE  L*u*v*  space  is  approximately  proportional  to  the 
Euclidean  distance  between  their  points.  Euclidean  distance  is  nothing  more 
than  our  common  sense  notion  of  distance  whose  formula  is  given  by: 

d(i,j)  +  (u|  -  uj}^  +  (vt  -  (1) 


wher  e 


cl(i,j)  =  the  distance  from  point  i  to  point  j 
(L|,  u^,  vt)  =  the  coordinates  of  point  i 
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Figure  2.  1976  CIE  UCS  Diagram 

Carter  and  Carter  (1982)  have  developed  a  computer  algorithm  to  solve  the 
color-spacing  problem  and  cite  several  military  applications.  They  describe 
how  the  method  could  be  used  to  choose  colors  for  strategic  aircraft  dis¬ 
plays  where  the  different  colors  correspond  to  friendly,  hostile,  and 
neutral  forces  or  various  enemy  target  types.  As  another  example,  an  air 
traffic  controller's  display  can  show  airplanes  at  different  altitudes,  all 
represented  by  various  colors.  Some  engineering  applications  include 
showing  "the  distribution  of  some  property  throughout  a  system,  such  as 
stress  in  a  structure  or  percent  of  full  capacity  in  various  parts  of  an 
electric  power  grid"  (Carter  and  Carter,  1982).  As  one  can  see,  there  are 
many  interesting  applications  of  the  color-spacing  problem,  both  in  the 
military  and  in  the  private  sector.  AFAMRL,  however,  has  not  been  able  to 
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implement  Carter  and  Carter's  program  and,  therefore,  has  not  used  optimiza¬ 
tion  techniques  when  decisions  regarding  color  selection  were  required. 

Concerning  the  algorithm.  Carter  and  Carter  themselves  believe  that  "presum¬ 
ably  a  more  efficient  one  could  be  devised  by  specialists  in  operations 
research"  (1982).  The  general  objectives  herein  are  to  generate  a  working 
program  that  produces  good  answers  to  the  color-spacing  problem,  and  is  also 
an  improvement  over  Carter  and  Carter's  method. 
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Section  3 

LITERATURE  REVIEW  OF  CARTER  AND  CARTER'S  METHOD 

The  only  known  solution  to  the  color-spacing  problem  was  Introduced  by 
Carter  and  Carter  (1982).  It  consists  of  first  randomly  placing  n  points  in 
the  region  (shown  in  Figure  1),  then  Identifying  the  minimum  CIE  L*u*v*  dis¬ 
tance  between  all  n(n-l)/2  pairs  of  points  In  that  region.  Let's  call  that 
value  D.  Once  that  distance  has  been  Identified,  the  two  closest  points 
(i.e.,  the  endpoints  associated  with  the  minimum  distance)  are  Investigated 
to  see  what  effect  is  created  by  moving  each  endpoint  to  Its  26  adjacent 
locations.  This  step  can  be  thought  of  as  "wiggling"  the  two  closest  points 
to  see  if  they  can  be  moved  farther  apart.  Figure  3  shows  the  various 
alternatives  for  one  endpoint.  The  alternatives  fall  on  the  boundaries  of  a 
cube  whose  sides  are  twice  the  step-size  In  length  and  whose  center  (the 
27th  point)  is  the  endpoint  Itself.  There  are  52  different  moves  that  can 
be  made  (26  for  each  endpoint).  If  a  move  decreases  the  distance  between 
those  two  points  or  pushes  one  of  them  outside  of  the  region,  that  alterna¬ 
tive  is  no  longer  considered.  The  remaining  alternatives  are  then  ranked 
according  to  how  much  they  Increase  the  distance  between  the  endpoints. 


Figure  3,  Alternative  Locations  for  an  Endpoint 
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If  the  highest  ranking  alternative  causes  an  increase  in  D,  that  move  is 
made.  Otherwise,  the  second  highest  ranking  alternative  is  considered,  and 
so  on,  until  one  of  the  alternatives  causes  an  increase  in  D.  If  none  of 
the  alternative  moves  increases  D,  then  the  step-size  is  halved  and  the  pro¬ 
cess  is  repeated;  otherwise,  the  point  is  moved  and  the  process  is  repeated 
until  an  expanding  move  cannot  be  made,  at  which  time  the  step-size  is 
hal ved. 

This  halving  continues  until  the  step-size  is  less  than  one  luminance  unit 
in  the  red,  green,  and  blue  color  space.  Once  the  step-size  is  less  than 
one  unit,  it  is  increased  back  to  its  original  value  and  the  point-moving 
and  step-size  halving  process  is  repeated,  "to  check  that  the  value  of  D 
arrived  at  is  not  merely  a  local  optimum"  (Carter  and  Carter,  1982).  If  the 
points  remain  in  the  same  place  throughout  the  repetition,  the  "configura¬ 
tion  of  points  is  assumed  to  represent  a  global  optimum  value  of  D"  (Carter 
and  Carter,  1982). 

In  their  results  section,  the  authors  state  that  when  the  number  of  colors 
(n)  is  three,  the  solution  is  identical  for  each  random  placement  of  points 
and,  in  general,  the  number  of  identical  solutions  decreases  as  n 
increases.  They  also  give  the  results  of  example  problems  where  3,  4,  6, 

10,  and  25  colors  were  spaced. 

Carter  and  Carter's  method  is  a  very  good  first  attempt.  It  also  has  some 
problems.  For  instance,  there  is  no  proof  that  the  solutions  are  globally 
optimal,  although  Carter  and  Carter  (1982)  assume  that  they  are.  In  the 
article,  the  authors  summarize  results  of  their  procedure  where  many  repli¬ 
cations  are  made  for  placement  of  3,  4,  6,  10  and  25  points.  For  the  case 
of  six  points,  50  replications  were  made  and  the  variance  among  the  values 
of  0  arrived  at  was  62.45,  where  the  maximum  value  of  0  for  all  replications 
was  124.08  (Carter  and  Carter,  1982).  This  relatively  high  variance  points 
to  the  fact  that  a  wide  range  of  values  for  D  can  be  expected  for  any  one 
replication,  and  that  most  replications  are  not  close  to  even  the  maximum 
known  value  of  D  for  that  problem,  let  alone  the  global  optimum  value  which, 
for  color-spacing  problems,  is  largely  still  unknown.  Whether  or  not  the 
solutions  represent  even  local  optima  is  uncertain  as  well. 
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Another  problem  is  that  there  is  no  validation  of  results.  For  instance. 
Carter  and  Carter  could  have  tested  the  algorithm  on  the  simple  case  of 
spacing  eight  po.nts  in  a  cube  to  see  if  the  optimal  solution  (a  point  in 
each  corner)  is  obtained.  There  is  not  even  a  discussion  as  to  the  physical 
desirability  of  the  colors  that  were  generated  by  the  algorithm. 

Also,  the  computer  code  is  poorly  documented  and  may  be  difficult  to  under¬ 
stand  since  variables  are  not  identified,  variable  names  are  often  conflict¬ 
ing  and  confusing,  and  statement  labels  are  poorly  numbered  and  sometimes 
not  found.  As  for  the  documentation,  the  number  of  internal  comments  is 
insufficient  to  make  the  program  easily  understandable. 

Carter  and  Carter's  method  has  its  good  points  as  well.  For  example,  it  is 
easy  to  understand.  Best  of  all,  the  technique  apparently  yields  solutions 
that  are  at  least  close  to  optimal,  which  is  all  they  really  set  out  to  do 
anyway.  But,  they  could  possibly  do  better.  For  instance,  the  two  closest 
points  may  be  optimally  spaced,  but  what  about  the  rest  of  them?  Wouldn't 
it  make  sense  to  go  through  the  algorithm  once,  then  maximize  the  second 
closest  distance,  the  third  closest,  and  so  on,  until  all  the  points  are 
optimally  spaced? 

Also,  Carter  and  Carter  maneuver  their  points  in  one  color  space  and  then 
transform  them  to  another  color  system  to  calculate  the  distances.  A  more 
appropriate  approach  would  be  to  do  all  the  maneuvering  and  calculations  in 
the  same  coordinate  system.  This  would  not  only  increase  the  efficiency  of 
the  program,  it  would  increase  its  flexibility  as  well.  It  would  enable  the 
:lgorithm  to  solve  any  type  of  spacing  problem,  not  just  the  color-spacing 
problem.  An  algorithm  like  Carter  and  Carter's,  but  with  better  validation 
and  documentation,  including  the  concept  of  successively  maximizing  the 
minimum  distance,  second  minimum  distance,  third  minimum,  and  so  on,  as  well 
as  continuous  operation  in  the  same  coordinate  system,  should  make  for  an 
improved,  more  reliable,  more  understandable,  more  efficient,  and  more 
flexible  solution  technique. 
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Section  4 

MATHEMATICAL  DEVELOPMENT 


The  first  part  of  this  section  will  show  the  mathematical  formulation  of  the 
color-spacing  problem  in  general  terms.  Next,  there  will  be  an  examination 
of  some  of  the  solutions  to  simple  spacing  problems.  The  final  segment  will 
discuss  how  the  problem  could  be  solved  by  nonlinear  programming. 

PROBLEM  FORMULATION 

For  the  spacing  problem,  we  want  to  maximize  the  minimum  distance  among  n 
points  in  a  convex  polyhedron.  A  polyhedron  can  be  simply  defined  as  a 
region  bounded  by  "many  faces."  A  triangle  is  an  example  of  a  polyhedron 
that  has  three  faces;  a  square  has  four  faces.  In  three  dimensions,  a  cube 
and  a  box  are  each  polyhedrons  with  six  faces  while  a  pyramid  has  five. 

These  are  all  convex  regions  as  well,  because  there  are  no  "juts"  sticking 
into  the  region.  A  star  (for  instance,  one  of  those  on  the  U.S.  flag)  is 
not  a  convex  region  for  that  very  reason.  More  formally,  a  region  is  convex 
if  every  point  in  that  region  can  be  connected  to  every  other  point  in  the 
region  with  a  straight  line  that  stays  completely  within  the  region.  All  it 
means  is  that  no  juts  or  dents  are  allowed. 

In  two  dimensions,  the  region  bounded  by  a  circle  is  convex,  as  is  that  of  a 
sphere  in  three  dimensions,  but  neither  is  a  polyhedron  because  they  have  no 
faces.  Mathematically  speaking,  a  polyhedron  must  have  linear  boundaries. 
That  means  that  a  polyhedron  must  be  made  up  of  a  series  of  straight  lines 
(or  planes  if  it's  three-dimensional).  No  curved  boundaries  are  allowed. 
Since  the  color-spacing  problem  is  three  dimensional,  we  will  primarily 
concern  ourselves  with  three  dimensions. 

The  boundaries  are  known  as  constraints.  A  linear  constraint  is  of  the 
form: 

ax  +  by  +  cz  =  k  (2) 
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where 


(x,y,z)  =  the  coordinates  of  a  point 
k  =  a  constant  known  as  the  right-hand-side  value 
a,  b,  and  c  =  any  real  numbers 

The  above  equation  describes  a  plane.  Changing  the  equation  to  ax  +  by  + 
cz  £  k,  we  have  the  points  on  a  plane,  in  addition  to  those  on  one  side  of 
it.  But  our  region  has  many  faces  to  it,  so  that  there  are  actually  many 
constraints  of  the  form: 

aix  +  bjy  +  cjz  _<  k^ 

a2X  +  bpy  +  C2Z  _<  k2 

aax  +  b3y  +  C3Z  _<  k3 

•  •  •  • 

•  •  •  « 

•  •  •  « 

a^x  +  b^y  +  CfZ  ^  k^r 

where 


f  =  the  number  of  faces  which  bound  the  region 

Together,  these  constraints  form  a  region  within  which  all  of  our  points 
must  stay.  If  a  point  is  outside  of  the  region,  then  at  least  one  of  these 
constraints  will  be  violated. 

Let's  turn  our  attention  again  to  the  objective.  We  want  to  maximize  the 
minimum  distance  among  n  points.  The  distance  from  one  point  to  another  in 
the  region  can  be  measured  by  the  Euclidean  distance  formula: 

=^(x,  -  +  (y,  -  yj)^  +  (z,  -  Zj)^  (3) 
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where 


cl(ij)  =  the  Euclidean  distance  from  point  i  to  point  j 
(xi,  y-j,  z-j)  =  the  coordinates  of  point  i 

In  order  to  solve  the  problem,  we  need  to  find  the  minimum  distance  between 
all  pairs  of  points.  If  there  are  only  three  points,  it's  easy;  there  are 
only  three  pairs  of  points  and  we  want  the  min  {d(l,2),  d(l,3),  d(2,3)} 
where  min  {  }  denotes  the  smallest  value  of  all  the  numbers  in  the  set  in 
brackets.  But  for  n  points,  that  set  may  be  quite  large.  The  minimum 
distance  for  n  points  would  be: 

min  {d(l,2),  d{l,3),  d(l,4),  ...,  d(l,n),  d{2,3),  d(2,4),  ..., 

d(2,n),  d(3,4),  ...,  d(3,n),  ...,  d(n-2,n-l),  d(n-2,n),  d(n-l,n)} 

As  a  matter  of  fact,  there  are  n(n-l)/2  values  in  the  set.  Letting 
D  =  min  (d(i,j)},  the  objective  then  is  to: 

maximize  D  =  maximize  min  {d(i,j)}  (4) 

subject  to 

=  d(1.2) 

=  d(1.3) 

=  d(n-l,  n) 
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^1^2  ^1^2  *  ‘-1^2  —  *^1 

•  •  •  • 

•  •  •  • 

•  •  •  • 

ai^n  +  +  CjZn  <  ki 

a2Xi  +  b2yi  +  C2Z1  ±  ^2 

®2^2  ^2'y2  ^  ^2^2  —  ^2 

•  •  •  • 

•  •  •  • 

•  •  •  « 

a2^n  +  ^  ^2^r.  1  '^2 


a^Xj  +  b^y^  +  c^Zj  ^  k^ 
afX2  +  bfy2  +  —  ^f 


^f’^n  +  ^fyn  ^  Cf^n  1  ''f 

where 

i  =  1  to  n-1 
j  =  i  +  1  to  n 
n  =  the  number  of  points 
f  =  the  number  of  faces 

the  first  set  of  constraints  =  the  distance  equations 
the  second  set  of  constraints  =  the  boundary  constraints  for 
all  n  points 
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To  uncomplicate  the  equations,  we  can  eliminate  the  square-root  signs  in  the 
distance  constraints  by  redefining  the  d(i,j)  to  be  distance  squared.  This 
will  not  affect  th“  fundamental  formulation.  Also,  we  can  subtract  D  from 
all  the  distance  constraints  to  create  equations  of  the  form; 

(Xi  -  Xj)^  +  (y^  ■  ^^i  -  Zj)^  -  D  =  d(i  ,j)  -  0  (5) 

But  we  know  that  d(i,j)  -  D  is  going  to  be  a  number  greater  than  or  equal  to 
zero  because  D  represents  the  smallest  possible  value  of  any  d(i,j).  There¬ 
fore,  we  can  change  the  equation  above  to: 

(x^.  -  Xj)2  +  (y.  -  y.]^  +  (z^  -  z.f  -  D  >  0  (6) 

The  problem  now  becomes: 

maximize  0  (7) 


subject  to 

(Xj  -  Xj)  (y,  -  yj)^  +  (z,  -  Zj]^  -  0  >  0 

Vg*  Vg*  Vgi''l 
>2*g  *  Vg  *  Vq  i  "2 


’f’-g  "  Vg  *  'f^g  i  "f 


where 


1  =  1  to  n  -  1 
j  =  i  +  1  to  n 

g  =  a  subscript  for  the  1  through  n  points  that  must  stay  within 
each  boundary  constraint 
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n  =  the  number  of  points 
f  =  the  number  of  faces 

There  are  a  few  more  simplifications  we  want  to  make.  For  practically  all 
commercial  methods  of  constrained  optimization,  the  constraints  are  required 
to  have  equal  signs  instead  of  inequalities.  We  accomplish  this  by  using 
variables  that  are  commonly  referred  to  as  surplus  and  slack  variables.  For 
example,  we  already  know  that  the  value  of  d(i,j)  -  D  is  not  negative,  but 
we  don't  know  how  large  it  is.  If  we  call  this  value  s^j  and  subtract  it 
from  the  distance  constraint,  we  have: 

(x,  -  Xj)2  *  (y,  -yj)2  t  (2,  -  z/  .  D  -  =  0  (8) 

The  constraint  is  now  an  equality  constraint  as  required,  and  s^j  is  called 
the  surplus  variable  because  it  represents  how  much  greater  than  zero 
equation  (6)  is. 

Similarly,  we  can  add  what  is  called  a  positive  slack  variable  to  the 
boundary  constraints  to  make  equations  of  the  form: 

%’'g  "  Vg  *  'h^g  *  Sg  =  "h  <’> 

where 

Shg  =  the  slack  variable 
The  problem  can  then  be  formulated  as: 

maximize  D  (10) 

subject  to 

(x,  -  Z.f  *  (y,  -  yj)^  t  (2,  -  2j)^  .  D  -  S,J  =  0 

’h-'g  *  Vg  *  S^g  *  ^hg  “  '‘h 
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where 


g’  ^g’  "g’ 


D,  s, 


ij’  hg’ 


i  -  1  to  n  -  1 
j  =  i  +  1  to  n 

g  =  a  subscript  for  the  1  through  n  points 

h  =  a  subscript  for  the  1  through  f  boundary  constraints 

n  =  the  number  of  points 
f  =  the  number  of  faces 

s^j  =  surplus  variables  for  the  distance  constraints 
S^g  =  slack  variables  for  the  boundary  constraints 

The  third  set  of  constraints  states  that  x,  y,  z,  D,  s,  S,  and  k  must  all  be 
greater  than  or  equal  to  zero,  commonly  referred  to  as  non-negativity 
constraints.  This  is  another  requirement  of  commercial  optimization  tech¬ 
niques.  We  already  know  that  D  is  positive  (negative  distance  is  impos¬ 
sible),  and  s^j  and  S^g  are  also  positive  as  discussed  before.  If  the  x's, 
y's,  and  z's  could  possibly  come  out  negative,  one  could  make  transforma¬ 
tions  of  the  form  x^  =  x^  -  t  where  t  is  the  minimum  possible  value  of  x,  to 

ensure  that  x^  is  positive.  If  k  is  negative,  both  sides  of  the  equation 
can  be  multiplied  by  -1  to  alleviate  the  problem,  but  the  sign  changes  from 
less  than  to  greater  than,  and  Sf,g  becomes  a  surplus  variable  instead  of  a 
slack  variable.  These  steps  ensure  that  S  and  k  will  always  be  positive 
and,  thus,  the  non-negativity  constraints  will  be  satisfied. 

Notice  that  we  now  have  n(n-l)/2  distance  constraints  that  are  nonlinear, 
and  n  •  f  linear  boundary  constraints.  The  next  step  is  to  incorporate  the 
nonlinear  constraints  into  the  objective  function  with  what  is  called  a 
Lagrangian  function.  If  we  call  the  left  side  of  the  first  distance  con¬ 
straint  Qj2,  the  second  0^3,  and  so  on,  the  Lagrangian  function  would  look 
like  this: 


D  +  XJ2Q12  ^13^13  ^ 


*  *lnOln* 


^  ^n-ln^n-ln 


(11) 
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or  similarly: 


D  + 


n-1 

E 

i=l 


n 

E 

j=i+l 


^ij 


[(»(  -  ♦  (y,  -  yjf  +  (2, 


where 


L  =  the  Lagrangian  function 
=  the  Lagrangian  variables 

So,  the  formulation  of  the  problem  now  looks  like  this: 


maxi mi ze 


n-1  n 

L  =  0>  L  E 

i=l  j=i+l 


ij 


[(x^  -  +  (y.  -  y.f  +  (z.  -  z.f 


subject  to 


3uX  +  b.  y  +  c.  z  +  S. 
h  g  h-'g  h  g  hg 


X,  y,  z,  D,  s,  S,  k  >  0 


where 


g  =  a  subscript  for  the  1  through  n  points 

h  =  a  subscript  for  the  1  through  f  boundary  constraints 

i  =  1  to  n  -  1 

j  =  i  +  1  to  n 

n  =  the  number  of  points 

f  =  the  number  of  faces 

s^j  =  surplus  variables  for  the  distance  constraints 
Sf,g  =  slack  variables  for  the  boundary  constraints 
(Xg,  yg,  Zg)  =  coordinates  of  point  g 
=  Lagrangian  variables 


We  now  have  a  nonlinear  objective  function  and  n  •  f  linear  constraints. 
Additionally,  the  objective  function  is  mathematically  classified  as  convex 
because  it  is  the  summation  of  many  convex  functions.  This  is  somewhat 
unfortunate  because,  if  it  were  concave,  the  problem  could  be  solved  using 
convex  programming;  a  solution  technique  which  is  widely  used  in  operations 
research.  However,  convex  programming  requires  convex  constraints  and  a 
concave  objective  function,  whereas  the  present  problem's  objective  function 
and  constraints  are  all  convex. 

Before  leaving  this  section,  one  more  thing  should  be  mentioned.  Origi¬ 
nally,  it  was  thought  that  the  color-spacing  problem  had  completely  linear 
boundary  constraints  (i.e.,  that  the  region  was  a  polyhedron).  It  turns  out 
that  this  is  not  the  case.  It  is  indeed  a  convex  polyhedron  in  the  red, 
green,  blue  coordinate  systems  but,  in  the  CIE  L*u*v*  system  (where  dis¬ 
tances  are  computed),  it  is  not.  The  reason  is  that  the  transformation  to 
CIE  L*u*v*  is  nonlinear  and  this  causes  the  faces  to  become  curved  sur¬ 
faces.  However,  the  region  is  still  convex  and  the  general  formulation 

above  still  applies.  Those  constraints  that  are  nonlinear  can  be  moved  into 
the  Lagrangian  function  in  much  the  same  manner  as  the  distance  constraints 
were,  while  the  linear  boundary  constraints  can  stay  put.  That  concludes 
the  problem  formulation,  and  we  are  now  ready  to  look  at  the  solutions  to 
some  easy  spacing  problems. 

SOLUTIONS  TO  SIMPLE  PROBLEMS 

Although  the  solutions  here  are  for  very  simple  problems  and  many  can  be 
recognized  intuitively  without  any  trouble,  some  are  not  so  apparent.  Shown 

here  are  the  solutions  for  spacing  2,  3,  4,  5,  and  6  points  in  a  square  to 

provide  a  feeling  for  the  types  of  solutions  that  can  be  expected  from  the 
spacing  problem.  Included  is  a  discussion  of  the  mathematical  aspects  of 
some  of  the  solutions  and  some  conclusions. 

Taking  a  look  first  at  the  solution  for  two  points  (shown  in  Figure  4),  note 
that  it  is  exactly  what  one  would  expect.  The  points  are  in  opposite  cor¬ 
ners.  Realize  also  that  there  are  actually  two  optimal  solutions.  One  is 
pictured  here  and  the  other  would,  of  course,  have  the  points  in  the  other 
two  corners. 
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Figure  4.  Solution  for  Spacing  Two  Points  in  a  Square 

The  solution  for  spacing  three  points  is  less  obvious  (see  Figure  5).  It 
turns  out  that  the  points  form  an  equilateral  triangle  (i.e.,  a  triangle 
whose  sides  are  of  equal  length)  with  one  point  in  a  corner  and  the  other 
two  points  on  opposite  faces  at  angles  of  15  degrees  from  the  corner 
point.  Not  coincidentally,  their  coordinates  are  (0,0),  (1,  tan  15  degrees) 
and  (tan  15  degrees,  1),  respectively,  for  a  square  whose  sides  are  one  unit 
long.  In  addition,  it  is  readily  apparent  that  this  problem  actually  has 
four  optimal  solutions;  each  one  corresponding  to  a  different  corner  point. 


Figure  5.  Solution  for  Spacing  Three  Points  in  a  Square 
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Spacing  four  and  five  points  in  a  square  is  trivial.  For  four,  the  points 
belong  in  the  corners.  The  five-point  solution  is  identical  to  the  four- 
point  solution  with  the  additional  point  going  in  the  center  of  the  region 
(see  Figures  6  and  7).  Incidentally,  both  of  these  problems  have  only  one 
optimal  solution. 


Figure  6.  Solution  for  Spacing  Four  Points  in  a  Square 


Figure  7.  Solution  for  Spacing  Five  Points  in  a  Square 

The  optimal  solution  for  spacing  six  points  in  a  square  is  shown  in 
Figure  8.  Most  people  would  probably  guess  that  all  four  corner  points 
would  be  filled  and  that  the  two  other  points  would  be  somewhere  in  the 
interior.  Instead,  there  are  two  corner  points,  one  interior  point,  and 
three  points  on  the  faces.  The  presence  of  a  counterintuitive  solution  like 
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Figure  8.  Solution  for  Spacing  Six  Points  in  a  Square 

this  one  emphasizes  the  need  for  a  computerized  solution  technique  like  the 
one  introduced  here.  Notice  that  the  left  two-thirds  of  the  region  is 
identical  to  a  skinny  five-point  solution,  and  that  the  points  in  the  right 
two-thirds  form  a  diamond  where  all  the  sides  are  of  equal  length.  In  all, 
there  are  six  distances  that  are  identical— the  four  of  the  diamond  and  the 
two  connecting  the  interior  point  to  the  corner  points--and  they  all  corre¬ 
spond  to  the  minimum  distance.  Notice  also  that  there  are  four  optimal 
solutions  for  this  problem. 

Let  us  consider  the  conclusions  that  can  be  drawn  from  these  simple 
examples.  First,  it  appears  that  the  points  generally  tend  to  fall  in  the 
corners  and  the  faces  before  they  go  to  the  interior  of  the  region.  Second, 
the  solutions  are  not  always  what  one  would  expect.  The  six-point  problem 
yields  a  good  example  of  that  phenomenon.  Third,  most  of  the  problems  have 
multiple  optimal  solutions,  so  that  the  decision  maker  can  often  choose 
between  several  alternatives.  Finally,  it  appears  that  the  optimal  solu¬ 
tions  all  have  a  high  percentage  of  their  distances  identical,  and  many  of 
these  are  equal  to  the  maximized  minimum  distance.  For  instance,  the  three- 
point  solution  has  3  out  of  the  3  equal  to  D,  the  four-point  solution  has  4 
out  of  6,  the  five-point  has  4  out  of  10,  and  the  six-point  has  6  out  of 
15.  This  means  that  one  might  be  able  to  tell  how  good  a  solution  is  not 
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only  by  the  value  of  D,  but  also  by  the  number  of  distances  that  are  equal 
to  that  value. 


SOLUTIONS  BY  NONLINEAR  PROGRAMMING 

Up  to  now,  you  may  have  been  thinking  that  the  problem  looks  interesting, 
but  the  mathematical  formulation  appears  complex.  How  does  one  actually 
solve  it?  There  are  several  ways.  One  may  write  a  heuristic  algorithm,  as 
Carter  and  Carter  did,  and  as  done  here.  Or,  one  may  solve  it  intuitively 
if  it's  simple  enough.  Or,  one  might  use  a  computerized  nonlinear  optimiza¬ 
tion  technique.  This  section  will  discuss  one  such  nonlinear  programming 
method. 

Originally  called  the  Method  of  Approximation  Programming  (MAP)  when  it  was 
introduced  in  the  late  1950s  by  Griffith  and  Stewart  (1961),  it  is  perhaps 
better  known  now  as  successive  linear  programming  (SLP).  It  works  like 
this:  Each  nonlinear  constraint  is  linearized  using  a  first-order  Taylor's 
series  expansion.  The  Taylor's  series  expansion  of  an  equation  is  a  mathe¬ 
matical  series  of  summed  terms  which  approximates  the  original  equation  and 
converges  to  the  original  as  terms  are  added  to  the  series.  Once  the  non¬ 
linear  constraints  are  linearized,  the  user  selects  a  starting  set  of  values 
and  solves  the  now  completely  linear  problem  using  linear  programming 
(LP).  The  solution  to  that  linear  programming  problem  is  then  used  as  the 
new  starting  set,  and  another  LP  solution  is  generated.  This  goes  on  until 
the  LP  solutions  become  identical,  at  which  time  the  LP  solution  is  also  the 
solution  to  the  nonlinear  problem. 

There  are  certain  restrictions  on  the  type  of  problem  for  which  this  tech¬ 
nique  will  guarantee  an  optimal  solution,  yet  Griffith  and  Stewart  (1961) 
state  that  "problems  have  been  solved  with  MAP  which  do  not  fully  satisfy 
all  of  these  requirements."  The  reason  this  technique  was  not  chosen  to 
solve  the  col  or- spacing  problem  is  because  of  the  size  of  the  problem.  For 
example,  in  order  to  solve  Carter  and  Carter's  problem  for  25  points,  there 
would  be  175  boundary  constraints  (some  linear,  some  nonlinear),  300  non¬ 
linear  distance  constraints,  and  76  decision  variables.  It  is  uncertain 
whether  an  LP  program  could  solve  that  problem  in  a  reasonable  amount  of 
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time,  let  alone  solve  it  again  for  the  second  closest  distance,  the  third 
closest,  and  so  on.  For  this  reason,  an  approach  using  a  heuristic  algo¬ 
rithm  was  chosen  instead  of  SLP. 


Section  5 

DESCRIPTION  OF  HEURISTIC 

The  heuristic  described  here  is  largely  patterned  after  that  of  Carter  and 
Carter,  which  was  described  in  Section  III.  However,  there  are  some  major 
differences.  The  first  portion  of  this  section  describes  how  the  algorithm 
works  and  some  of  the  logic  behind  it.  The  second  portion  identifies  the 
differences  between  this  algorithm  and  Carter  and  Carter's. 

THE  ALGORITHM 

Recall  that  the  purpose  of  this  effort  was  to  maximize  the  minimum  distance 
among  n  points  in  a  convex  region,  and  then  to  successively  maximize  the 
second  minimum  distance,  the  third  minimum,  and  so  on,  until  all  points  are 
optimally  spaced.  The  first  step  in  achieving  this  objective  is  to  randomly 
place  the  points  uniformly  throughout  the  region.  For  the  color-spacing 
problem,  the  region  looks  like  the  one  shown  in  Figure  9,  Recall,  from 
page  22,  that  this  region  is  not  really  a  polyhedron,  although  Figure  9 
depicts  it  as  one.  The  actual  color-spacing  region  would  have  many  of  its 
boundaries  bowed  out  slightly.  Because  they  are  bowed  out  and  not  in,  the 
region  is  still  convex,  so  the  solution  technique  described  here  still 
appl i es. 

In  step  2,  the  program  finds  the  minimum  distance  among  all  n(n-l)/2  pairs 
of  points,  and  the  endpoints  that  correspond  to  that  distance.  The  third 
step  is  to  define  the  26  alternative  locations  for  each  of  the  two  endpoints 
found  in  step  2.  These  alternatives  lie  on  a  unit  cube  around  the  endpoint 
whose  sides  are  twice  the  step-size  in  length.  Figure  3  (page  11)  depicts 
these  alternative  endpoint  locations. 

In  the  fourth  step,  the  alternatives  are  checked  to  see  if  an  improvement 
can  be  made  by  moving  the  endpoints.  The  program  does  this  by  first  check¬ 
ing  the  alternatives  to  see  if  they  are  in  the  region.  Then,  taking  one 
legal  alternative  at  a  time,  it  checks  the  distance  between  endpoints  to  see 
if  that  value  is  greater  than  before.  Next,  it  checks  the  new  overall  mini¬ 
mum  distance  to  see  if  that  is  greater.  If  all  of  these  conditions  are  met. 
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Figure  9.  Approximate  Feasible  Region  for  the  Color  Spacing 
Problem  in  the  CIE  L*u*v*  Color  System 

the  endpoints  are  moved  and  the  routine  checks  the  next  alternative  for  an 
even  greater  possible  improvement.  It  does  this  for  every  possible  combi¬ 
nation  of  alternative  endpoint  locations;  a  total  of  729  combinations.  If 
several  alternatives  have  the  same  overall  minimum  distance,  which  is 
likely,  the  tie  goes  to  the  alternative  producing  the  greatest  distance 
between  endpoints. 

The  program  then  returns  to  step  2  and  finds  the  new  overall  minimum  dis¬ 
tance  and  the  endpoints  that  correspond  to  that  distance.  This  loop  con¬ 
tinues  until  no  improvement  can  be  made  with  the  given  step- size.  At  that 
time,  the  step-size  is  halved  and  the  program  begins  looping  again,  starting 
at  step  2.  This  halving  process  continues  until  the  step- size  reaches  some 
minimum  value  chosen  by  the  user  and  stej?  5  is  started. 

The  fifth  step  controls  the  fixing  of  points  already  maximized  and  the 
status  of  the  program.  When  the  step-size  reaches  its  minimum  tolerance, 
one  of  the  endpoints  is  fixed  in  its  position,  the  step-size  is  returned  to 
its  original  value,  and  the  program  returns  to  step  2.  The  only  difference 
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in  the  logic  now  (other  than  a  point  being  fixed)  is  that  the  minimum  dis¬ 
tance  calculated  in  step  2  can't  have  two  fixed  points  as  its  endpoints. 

This  enables  the  program  to  concentrate  on  the  second  shortest  distance. 

The  fixing  of  points  continues  until  all  the  points  become  fixed,  at  which 
time  they  are  all  unfixed  and  the  routine  starts  all  over  again  at  step  2  as 
if  the  current  location  scheme  were  the  starting  positions.  If  the  program 
runs  all  the  way  through  again  and  the  points  haven't  moved,  then  the  pro¬ 
gram  is  done.  Otherwise,  the  program  keeps  going  until  there  is  no  change 
in  location  for  one  complete  iteration.  Figure  10  contains  a  flowchart  of 
the  program's  logic,  while  Appendix  C  contains  an  example  of  how  the  program 
optimally  spaced  four  points  in  a  square. 

DIFFERENCES  WITH  CARTER  AND  CARTER'S  METHOD 

As  stated  at  the  beginning  of  this  section,  there  is  a  resemblance  between 
Carter  and  Carter's  method  and  the  one  described  here.  However,  there  are 
four  major  differences. 

The  first  of  these  differences  lies  in  the  checking  of  alternatives  to  see 
what  expanding  move  should  be  made.  Carter  and  Carter  move  only  one  end¬ 
point  at  a  time  and,  therefore,  choose  between  only  52  possible  endpoint 
location  schemes  (26  for  each  endpoint).  On  the  other  hand,  the  Roley  algo¬ 
rithm  moves  both  endpoints  at  once  and,  as  a  result,  has  a  total  of  729  pos¬ 
sible  endpoint  location  schemes.  This  constitutes  a  more  exhaustive  search 
for  the  best  possible  move. 

Along  the  same  lines.  Carter  and  Carter's  method  ranks  the  alternatives  as 
to  how  much  they  increase  the  distance  between  endpoints  and  then  chooses 
the  highest  ranking  alternative  that  also  increases  the  overall  minimum 
distance.  On  the  other  hand,  the  Roley  program  chooses  the  alternative  that 
most  increases  the  overall  minimum  distance,  with  ties  going  to  the  alterna¬ 
tive  which  causes  the  greatest  increase  in  the  endpoint  distance.  The 
difference  is  that  Carter  and  Carter's  algorithm  may  not  choose  the  best 
alternative  for  increasing  the  overall  minimum  distance,  whereas  the  Roley 
program  does. 
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Figure  10.  Flowchart  of  Roley  Heuristic 


Another  difference  has  to  do  with  the  thoroughness  of  the  spacing  method. 
Carter  and  Carter  maximize  the  minimum  distance  but  leave  the  second  mini¬ 
mum,  third  minimum,  and  all  the  others  where  they  are.  The  Roley  algorithm 
not  only  maximizes  the  minimum  distance,  it  successively  maximizes  the  sec¬ 
ond  minimum  distance,  the  third  minimum,  and  so  on  until  all  points  are 
optimally  spaced.  The  result  is  a  complete  optimization  of  all  distances, 
not  just  the  minimum  distance.  The  price  is  a  more  complex  and  slower  com¬ 
puter  algorithm. 

The  fourth  and  final  difference  has  to  do  with  the  respective  coordinate 
systems  that  are  used  as  the  primary  operating  spaces.  Carter  and  Carter's 
technique  spaces  points  in  the  red,  green,  and  blue  color  system,  and  con¬ 
verts  the  colors  to  CIE  L*u*v*  color  coordinates  i'  jrder  to  calculate 
spatial  distances.  The  Roley  program  acts  in  a  reverse  manner.  It  maneu¬ 
vers  points  and  calculates  distances  in  the  CIE  L*u*v*  system,  but  checks 
color  locations  in  the  red,  green,  and  blue  system  to  make  sure  they  are 
within  the  boundaries  of  the  region.  It  was  originally  hoped  that  all 
operations  and  calculations  could  be  done  in  the  CIE  L*u*v*  system  for 
increased  efficiency,  but  the  equations  for  the  boundary  constraints  could 
not  be  derived  in  the  L*u*v*  system,  preventing  its  use  as  a  location¬ 
checking  coordinate  system. 
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Section  6 
RESULTS 


There  are  several  questions  that  need  to  be  answered  in  this  sectiori.  Oce: 
the  algorithm  guarantor*  an  opti'^al  solution?  Does  it  work  for  very  Ic.-ge 
problems?  How  well  does  it  perform  in  comparison  to  Carter  and  Carter's 
method?  What  are  some  of  the  factors  to  which  the  program  is  sensitive? 

What  are  its  biggest  advantages?  What  are  its  biggest  disadvantages? 

optimality 

The  algorithm  does  not  guarantee  an  optimal  solution.  However,  it  was 
tested  on  a  number  of  problems  whose  solutions  are  known  and  the  results  are 
encouraging.  These  problems  include  spacing  3,  4,  5,  and  6  points  in  a 
square  and  spacing  eight  points  in  a  cube.  Each  problem  was  tested  with 
three  different  sets  of  randomly  placed  points  and  the  algorithm  produced 
the  optimal  solution  in  14  of  the  15  problems.  That's  a  success  rate  of 
93  percent.  The  only  unsuccessful  attempt  involved  one  of  the  tries  at 
spacing  five  points  in  a  square.  In  that  instance,  the  minimum  distance 
arrived  at  was  .638  units  compared  with  the  optimal  solution  of  .707,  so 
that  it  was  within  90  percent  of  optimal.  In  general,  the  more  replications 
that  are  performed,  the  better  the  chance  of  success. 

SIZE  LIMITATIONS 

This  may  be  the  largest  downfall  of  the  Roley  program.  The  original  inten¬ 
tion  was  for  it  to  accommodate  up  to  50  points  and  10  faces.  However,  the 
largest  problem  it  has  been  able  to  successfully  solve  was  the  color-spacing 
problem  with  23  colors  and  7  faces.  Larger  problems  exceeded  the  time  limit 
of  1000  CPU  seconds  on  the  CDC  845  (Cyber)  computer.  The  reason  for  this 
size  deficiency  boils  down  to  the  thoroughness  of  the  program.  Each  itera¬ 
tion  checks  n(n-l)/2  distances  for  up  to  729  alternatives.  It  takes  21 
iterations  for  the  step-size  to  reach  its  minimum  value.  The  program  must 
perform  those  21  iterations  for  each  of  the  n  points  that  are  fixed  in  suc¬ 
cession.  It  then  goes  through  the  whole  process  at  least  one  more  time  or 
until  the  point  locations  do  not  change.  Experience  has  shown  that  it 
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usually  goes  through  twice.  With  an  increase  in  n,  the  number  of  computa¬ 
tions  (and  thus  the  CPU  time  required)  goes  up  exponentially,  so  one  can 
easily  tell  that  size  is  a  serious  factor.  Presumably,  reducing  the  number 
of  alternatives,  increasing  the  minimum  step-size,  or  reducing  the  amount  of 
point  fixing  would  increase  the  ability  of  the  program  to  handle  larger 
problems. 

PERFORMANCE  IN  COMPARISON  TO  CARTER  AND  CARTER'S  METHOD 

The  comparison  of  this  program  with  Carter  and  Carter's  was  done  using 
identical  for  Y^,  Uq,  Vq,  and  assuming  an  identical  CRT.  Carter  and  Carter 
tested  their  program  using  a  Univac  1182  computer  running  a  batch  operating 
system.  Originally,  the  Roley  program  was  tested  on  a  CDC  845  computer, 
also  running  a  batch  operating  system.  Shortly  before  the  present  report 
was  published,  however,  additional  results  became  available  for  the  Roley 
program,  running  on  a  DEC  POP-11/60  minicomputer.  The  PDP-ll  in  question 
has  an  FP-11  floating-point  processor,  runs  the  RSX-ll/M  timesharing  oper¬ 
ating  system  (version  4.1,  patch  D),  and  has  a  FORTRAN-77  compiler  (RSX 
version  5.0),  which  was  used  to  compile  the  program.  Table  1  contains  com¬ 
parative  results  for  spacing  3,  4,  6,  10,  and  25  colors.  The  times  shown 
for  the  PDP  are  not  pure  execution  times  but  are  the  differences  between  the 
system  clock  times  at  start  versus  end-of- replication.  Thus,  they  reflect 
some  system  overhead  and  time  slicing.  During  program  execution,  other 
tasks  on  the  system  consisted  entirely  of  editing  and  occasional  text  for¬ 
matting  runs.  Therefore,  the  times  are  representative  of  what  should  be 
expected  for  a  lightly  loaded  CPU. 

The  original  intention  was  to  perform  only  five  replications  for  each  prob¬ 
lem  because  time  considerations  prevented  performing  50  runs,  as  Carter  and 
Carter  did.  However,  20  runs  of  the  ten  color  problem  were  made  on  the  CDC 
in  an  attempt  to  decrease  the  disparity  of  results.  The  25-color  problem 
could  not  be  executed  on  the  CDC  because  the  CPU  limit  of  1000  seconds  was 
exceeded  on  each  attempt.  Because  the  execution  time  on  the  PDP  for  the 
25-color  problem  was  so  excessive,  only  one  replication  was  performed.  How¬ 
ever,  as  can  be  seen,  the  resulting  minimum  distance  for  that  replication  is 
slightly  better  than  Carter  and  Carter's. 
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TABLE  1.  COMPARISON  OF  MINIMUM  DISTANCES 


Number 

of 

Colors 

Method 

Best  Answer 
(CIE  L*u*v*) 

Variance 

Number 
of  Repli¬ 
cations 

Average 
Computer  Time 
Per  Replica¬ 
tion  (seconds) 

3 

Roley -CDC 

239.26 

1.17 

5 

23.50 

Roley-POP 

239.25 

395.76 

5 

66.60 

Carter-Uni  vac 

239.33 

.11 

49 

1.31 

4 

Roley-CDC 

156.28 

136.62 

5 

42.95 

Roley -PDP 

156.10 

119.66 

5 

156.00 

Carter-Uni  vac 

155.87 

37.14 

50 

1.76 

r- 

D 

Roley-CDC 

124.85 

79.61 

5 

100.49 

Roley-PDP 

124.47 

9.34 

5 

288.00 

Carter-Uni  vac 

124.08 

62.45 

50 

2.48 

10 

Roley-COC 

71.07 

48.61 

20 

178.90 

Roley-PDP 

90.01 

8.15 

5 

1628.60 

Carter-Uni  vac 

89.43 

27.73 

50 

4.22 

25 

Roley-CDC 

■■ 

••• 

5 

>1000.00 

Roley-PDP 

54.42 

— 

1 

-2  days 

Carter-Univac 

51.60 

7.34 

46 

12.39 

It  is  difficult  to  compare  this  algorithm  with  Carter  and  Carter's  because 
there  are  so  many  factors  that  can  be  used  as  measures  of  effectiveness. 
Most  people,  however,  would  consider  the  minimum  distances  produced  by  each 
method  to  be  of  primary  importance.  Notice  in  Table  1  that  most  of  the 
minimum  distances  are  similar,  although  Carter  and  Carter's  results  are 
based  on  many  more  runs.  The  minimum  distance  for  spacing  10  colors  on  the 
COC  demonstrates  the  importance  of  performing  a  sufficient  number  of  repli¬ 
cations.  Originally,  this  result  was  interpreted  as  indicating  a  possible 
deficiency  in  the  Roley  program  for  problems  involving  more  than  six 
points.  However,  when  the  POP  results  became  available.  It  became  clear 
that  this  failure  was  simply  a  matter  of  bad  luck  in  the  initial  random 
placement  of  points.  As  for  the  differences  in  the  variance  of  the  Roley 
program's  solutions  on  the  COC  versus  POP,  it  is  uncertain  whether  they 
reflect  a  true  difference  or  merely  the  small  sample  sizes. 
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Another  criterion  that  might  be  used  to  judge  the  two  methods  would  De  to 
compare  the  distances  between  the  second  closest  points,  third  closest,  and 
so  on.  This  would,  indeed,  be  an  appropriate  measure  because  the  funda¬ 
mental  conceptual  ditrerence  between  the  two  methods  is  that  the  Roley  algo¬ 
rithm  concentrates  on  successively  maximizing  all  distances,  whereas  Carter 
and  Carter's  maximizes  only  the  minimum  distance.  Presumably,  these  efforts 
should  have  paid  off,  but  it  is  difficult  to  judge  because  of  a  limited 
sample  size.  Table  2  presents  a  comparison  of  distances  for  the  six-color 
spacing  problem.  Although  it  would  be  desirable  to  compare  other  problems 
in  addition  to  the  six-color  one.  it  is  the  only  one  published  and  avail¬ 
able.  Notice  that  the  first  three  distances  are  slightly  better  for  the 
Roley  algorithm,  the  next  six  are  fairly  even,  and  the  last  six  are  clearly 
in  the  Carters'  favor.  This  could  indicate  that,  while  maximizing  the 
second  and  third  minimum  distances  produces  better  results  for  those  values, 
it  may  cost,  in  the  long  run,  in  the  form  of  smaller  values  for  the  greater 
distances.  Notice  also  that  the  first  through  sixth  distances  are  basically 
the  same  regardless  of  whether  or  not  they  are  maximized.  They  seem  to  have 
been  automatically  maximized  by  maximizing  the  minimum  distance.  This  indi¬ 
cates  that  successively  maximizing  greater  distances  may  not  be  worth  the 
extra  computer  time  required.  These  results  also  support  the  observation 
that  was  made  in  Section  IV,  i.e.,  that  optimal  solutions  will  generally 
have  several  distances  equal  to  the  minimum  distance.  This  is  a  charac¬ 
teristic  that  could  be  extremely  useful  in  determining  whether  or  not  a 
particular  solution  is  close  to  optimal. 

TABLE  2.  COMPARISON  OF  ALL  DISTANCES  FOR  THE  SIX-COLOR  PROBLEM 


ith 

Minimum 

Di stance 

Rol  ey 
CDC 

ith 

Mi nimum 
Carter  Distance 

Roley 

CDC 

Carter 

ith 

Minimum 

Di stance 

Roley 

CDC 

Carter 

1st 

125* 

124 

6th 

129 

130* 

nth 

182 

218* 

2nd 

125* 

124 

7th 

136* 

133 

12th 

201 

237* 

3rd 

125* 

124 

8th 

138 

140* 

13th 

220 

239* 

4th 

125 

125 

9th 

149* 

147 

14th 

239 

249* 

5th 

125 

125 

10th 

171 

210* 

15th 

241 

263* 

*Denotes  the  superior  value. 
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One  might  also  be  impressed  with  the  Roley  method  if  it  could  take  Carter 
and  Carter's  best  solution  for  a  particular  problem,  use  it  as  a  starting 
location  scheme,  and  improve  upon  it.  Using  Carter  and  Carter's  six-color 
solution  (because  it  was  the  only  one  available)  as  an  initial  set  of  colors 
resulted  in  a  mild  improvement  of  the  minimum  distance  from  124.026  CIE 
L*u*v*  units  to  124.145  units  on  the  CDC.  Research  by  Carter  and  Carter 
(1982)  has  shown  that  distinguishability  between  colors  "deteriorates 
rapidly  when  the  distance  between  colors  is  about  40  CIE  L*u*v*  units,"  so 
this  slight  improvement  of  .12  unit  is  inconsequential.  Perhaps  more  impor¬ 
tantly,  my  algorithm  improved  ai  11  out  of  the  14  other  distances  with  one 
unchanged  and  two  decreasing  slightly.  Note  that  the  original  distance  of 
124.026  does  not  correspond  to  Carter  and  Carter's  calculated  value  of 
124.08  as  shown  in  Table  1.  Taking  their  same  set  of  color  locations  in 
red,  green,  and  blue  coordinates,  a  minimum  distance  of  124.026  in  L*u*v* 
coordinates  instead  of  124.08  (the  exact  L*u*v*  coordinates  were  not  pub¬ 
lished)  was  obtained.  The  difference  is  probably  due  to  roundoff  error. 

The  improvement  on  their  optimal  solution  with  basically  the  same  teennique 
is  due  to  differences  in  accuracy.  Because  Carter  and  Carter's  minimum 
step-size  is  one  unit,  they  are  not  able  to  attain  the  same  level  of  accu¬ 
racy  as  the  Roley  program  does  with  a  minimum  step-size  of  .0001  units. 
Consequently,  the  Carters  are  not  able  to  "snuggle"  their  points  into  the 
best  locations.  All  this  indicates  is  that  the  minimum  step-size  is  a 
control  over  the  level  of  accuracy  one  wishes  to  achieve. 

As  was  briefly  discussed  before,  the  question  of  computer  time  may  also  be 
an  important  performance  characteristic  by  which  to  judge  the  two  methods. 
Theoretically,  Carter  and  Carter's  method  should  be  considerably  faster 
because  it  maximizes  only  the  minimum  distance  and  isn't  concerned  with  the 
second  minimum  distance,,  third  minimum,  and  so  on.  Looking  again  at 
Table  1,  one  can  see  that  this,  in  fact,  appears  to  be  the  case.  It  should 
be  remembered,  however,  that  comparing  ryn  times  on  different  computers  is 
like  comparing  apples  and  oranges.  The  CDC  versus  POP  results  show  this 
clearly. 

In  the  present  author's  opinion,  the  chief  criterion  for  comparison  is  which 
program  works.  The  present  author  knows  of  no  place  where  Carter  and 
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Carter's  program  is  currently  running  (other  than  their  own  facility)  and 
several  places  where  people  have  tried  to  implement  it  and  have  failed.  On 
the  other  hand,  the  Roley  program  is  two  for  two  so  far,  so  it  must  be  given 
the  edge.  Interested  parties  will  be  pleased  to  find  that  it  is  heavily 
documented,  is  easy  to  understand,  utilizes  classic  structured  programming 
concepts,  and  is  written  in  FORTRAN  77.  These  characteristics  give  it  the 
added  advantage  of  being  easier  to  modify  and/or  upgrade  by  others  who  may 
be  so  inclined.  A  listing  of  the  code  is  contained  in  Appendix  B  for  all  to 
judge  for  themselves. 

Briefly  summarizing  this  section,  we  found  that  Carter  and  Carter's  method 
works  equally  well  for  finding  the  maximum  minimum  distance,  but  requires 
more  runs.  For  the  six-color  spacing  problem,  theirs  yielded  slightly  worse 
values  for  the  second  and  third  minima,  and  better  results  for  the  greater 
distances.  Carter  and  Carter's  program  took  less  time  on  a  Uni  vac  1182,  but 
used  a  larger  minimum  step-size  and  produced  a  slightly  lower  degree  of 
accuracy.  Unlike  Carter  and  Carter's,  the  Roley  program  works  for  other 
users,  is  well  documented,  and  easy  to  understand. 

SENSITIVITY  TO  PARAMETERS 

The  program  introduced  in  this  report  is  sensitive  to  four  parameters.  They 
are  the  maximum  step-size,  the  step-size  reducing  scheme,  the  minimum  step- 
size,  and  the  initial  random  positioning.  The  maximum  step-size  is  impor¬ 
tant  because  some  solutions  are  not  obtainable  for  certain  step-sizes  and 
point  locations.  Take,  for  instance,  the  case  of  spacing  eight  points  in  a 
cube  whose  sides  are  one  unit  long.  Suppose  seven  of  the  points  are  in  the 
corners  and  the  eighth  is  in  the  very  middle  of  the  cube.  The  eighth  point 
should  migrate  to  the  empty  corner,  but  unless  the  step-size  is  greater  than 
one-third  of  a  unit,  the  minimum  distance  will  decrease  if  it  tries  to  move 
in  any  direction.  Essentially,  the  point  is  "locked  in"  to  its  location. 

To  avoid  this  problem,  the  use  of  a  maximum  step-size  equal  to  approximately 
half  the  distance  between  the  two  farthest  points  in  the  region  is  recom¬ 
mended.  This  way,  a  point  located  in  the  middle  can  move  to  any  other  spot 
in  the  region. 
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A  point  doesn't  have  to  be  in  the  middle  to  be  locked  in.  Suppose  thit  the 
eighth  point  in  our  example  above  is  not  located  in  the  exact  midJl  .  Sup¬ 
pose,  instead,  that  it  is  located  near  the  empty  corner,  but  less  chan  one- 
third  of  a  unit  from  the  middle.  A  step-size  of  .5  now  puts  the  point  out 
of  the  region,  so  it  still  can't  move  closer  to  where  it's  supposed  to  be. 

If  the  step-size  is  halved,  the  point  still  may  not  be  able  to  go  toward  the 
corner.  But,  some  other  step-size  reduction  scheme,  like  reducing  it  by 
four-fifths,  might  send  the  point  to  the  optimal  solution.  In  that  sense, 
the  reduction  scheme  is  another  factor  which  could  affect  the  solution.  The 
problem  with  that  alternative  is  that  it  takes  quite  a  few  more  iterations. 
For  instance,  to  go  from  .5  to  .0001  by  halves  takes  13  iterations.  To  go 
from  .5  to  .0001  by  four-fifths  takes  39  iterations.  More  iterations  means 
more  computation  time  and,  thus,  a  less  desirable  program.  Both  the  Roley 
and  the  Carter  and  Carter  programs  use  a  halving  scheme,  but  it  is  recom¬ 
mended  that  users  experiment  with  this  parameter  to  find  out  what  is  best 
for  any  given  application. 

The  minimum  step-size,  although  a  factor  to  the  solution,  is  not  as  impor¬ 
tant  as  the  two  parameters  discussed  above.  The  minimum  step-size  controls 
the  accuracy  of  the  solution.  For  three  decimal  places  of  accuracy,  .0001 
should  be  the  minimum  step-size.  For  two  places,  use  .001,  and  so  on. 
Increasing  the  minimum  step-size  can  cause  a  surprising  decrease  in  the  num¬ 
ber  of  iterations  required.  For  instance,  it  takes  only  six  iterations  to 
go  from  .5  to  .01  by  halves  as  opposed  to  13  iterations  to  go  to  .0001. 

The  final  parameter  that  significantly  affects  the  solution  is  the  seed. 

This  is  the  starting  value  that  is  used  in  the  random  number  generating 
sequence  which  calculates  the  starting  locations  for  the  points.  One  can 
intuitively  realize  the  importance  of  the  seed  by  understanding  that  some 
initial  location  sets  are  better  than  others.  A  bad  location  set  is  shown 
in  Figure  11.  None  of  the  points  in  this  example  can  move  anywhere  without 
getting  closer  to  another  point,  and  yet  the  solution  is  not  optimal.  It  is 
called  a  "locally"  optimal  solution  as  opposed  to  the  "globally"  optimal 
solution  which  has  a  point  in  each  corner.  Spacing  problems  apparently  tend 
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Figure  11.  Example  of  Locally  Optimal  Solution 

to  have  many  local  optima.  Changing  the  seed  and  doing  numerous  replica¬ 
tions  will  increase  the  chances  of  finding  a  global  optimum. 


ADVANTAGES 


The  heuristic  algorithm  described  here  has  numerous  advantages.  Some  of 
them  are: 


1.  Flexibility.  The  program  works  for  any  spacing  problem,  including 
the  color-spacing  problem.  It  has  a  myriad  of  uses  in  the  field 
of  optics,  as  described  by  Carter  and  Carter  (2:2939).  It  can 
even  be  extended  to  four  dimensions  for  possible  use  in  the  field 
of  physics.  It  can  accommodate  distance  formulas  other  than  the 
Euclidean  distance  formula  or  convex  regions  besides  poly¬ 
hedrons.  Its  uses  are  limited  only  by  the  user's  imagination. 

2.  Simplicity.  It  is  simple  to  understand  and  simple  to  use.  The 
structured  FORTRAN  programming  format  makes  it  simple  to  modify  or 
debug. 
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3.  Reliability.  It  has  been  tested  favorably  against  simple  problems 
with  known  solutions  and  more  complex  problems  like  Carter  and 
Carter's  color-spacing  problems.  It  has  proven  to  produce  good 
answers  with  fewer  replications  than  previous  solution  techniques. 

4.  Successive  Maximization.  It  not  only  maximizes  the  minimum 
distance,  it  successively  maximizes  the  second  minimum,  third 
minimum,  and  so  on.  Because  all  points  and  their  distances  in 
relation  to  each  other  are  important  in  the  spacing  problem,  this 
is  a  more  appropriate  objective. 

DISADVANTAGES 

Here  is  a  list  of  the  method's  major  disadvantages: 

1.  Optimality.  The  method  does  not  guarantee  a  globally  optimal 
solution,  nor  does  it  guarantee  even  a  locally  optimal  solution. 

It  just  guarantees  a  "good"  answer  and  it  may  take  several  runs  to 
get  that  "good"  answer. 

2.  Size  Limitations.  The  program  has  been  proven  on  a  problem  of  at 
most  25  points  and  seven  faces.  There  is  no  guarantee  that  it 
will  work  on  very  large  problems.  Computer  time  is  the  key 
resource  here  because  it  increases  exponentially  as  the  number  of 
points  increase.  It's  up  to  the  user  to  weigh  the  value  of  com¬ 
puter  time  against  the  value  of  a  solution. 
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Section  7 
SUMMARY 


Let  us  review  what  has  been  discussed  thus  far.  After  being  introduced  to 
the  spacing  problem,  we  learned  of  some  examples  in  the  real  world,  most 
notably  the  color-spacing  problem.  We  were  educated  on  the  theory  behind 
the  color-spacing  problem  and  we  found  out  how  Carter  and  Carter  proposed  to 
solve  it.  Next,  we  learned  how  the  spacing  problem  is  formulated  mathemati¬ 
cally  and  were  able  to  see  some  solutions  to  simple  problems.  We  also  found 
out  how  the  problem  could  be  solved  using  successive  linear  programming. 
Next,  we  were  introduced  to  a  new  heuristic  algorithm,  very  similar  to 
Carter  and  Carter's,  that  is  designed  to  solve  the  spacing  problem. 

Finally,  we  found  out  how  well  the  new  algorithm  works.  It  is  now  time  to 
identify  our  conclusions  and  some  suggestions  for  further  research. 

CONCLUSIONS 

The  program  works.  That  was  the  main  objective  of  this  effort  and,  obvi¬ 
ously,  it  was  accomplished.  Perhaps  equal  to  that  was  the  objective  of 
providing  AFAMRL  with  a  computerized  solution  to  the  color-spacing  prob¬ 
lem.  Another  objective  was  to  improve  upon  Carter  and  Carter's  method.  It 
is  difficult  to  compare  the  methods  because  time  constraints  prevented  doing 
a  full  statistical  analysis.  But,  if  one  considers  that  the  new  program 
runs,  is  easy  to  understand  and  use,  and  yields  solutions  that  are  compar¬ 
able  to  Carter  and  Carter's,  one  must  consider  it  an  improvement. 

The  present  author  originally  thought  that  Carter  and  Carter's  algorithm  was 
a  rough  first  attempt  generating  suboptimal  solutions  which  could  easily  be 
improved  upon.  It  is  a  credit  to  their  work  that  those  solutions  are  appar¬ 
ently  not  so  suboptimal.  But  it  must  be  remembered  that,  regardless  of 
which  method  is  used,  both  require  numerous  replications  to  come  up  with  an 
answer  that  is  only  guaranteed  to  be  "good,"  So  there  is  still  plenty  of 
room  for  improvement. 
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SUGGESTIONS  FOR  FURTHER  RESEARCH 


There  are  two  directions  that  subsequent  research  can  take.  One  is  to 
refine  the  new  algorithm.  The  other,  and  more  exciting  alternative,  is  to 
solve  the  spacing  problem  using  some  form  of  nonlinear  programming.  The 
method  described  in  Section  IV  of  this  report  has  some  definite  possibili- 
tips,  b'lt  the  difficulty  is  that  the  constraints  of  the  color  region  are  not 
linear  and  are  difficult  to  derive  mathematically.  Carter  and  Carter  (1982) 
describe  the  region  as  "approximately  a  triangular  prism,"  but  nowhere  in 
the  literature  is  there  a  mathematical  formulation  of  the  boundaries.  This 
could  prove  to  be  a  major  stumbling  block  for  anyone  attempting  to  solve  the 
color-spacing  problem  using  nonlinear  programming. 

Alternatively,  there  are  improvements  that  could  be  made  to  the  Roley  algo¬ 
rithm.  The  most  worthwhile  appears  to  be  the  use  of  vectors  and  matrices  to 
help  define  possible  new  positions  for  the  two  closest  points.  The  present 
method  disregards  any  alternative  that  takes  a  point  out  of  the  feasible 
region.  A  method  could  be  devised  that  would,  instead,  take  the  point  to 
the  boundary  of  the  feasible  region  and  still  consider  it  as  a  viable  alter¬ 
native.  Indeed,  it  may  be  a  desirable  move  because  of  the  tendency  of  the 
points  to  migrate  to  the  boundaries  and  corners  of  the  region.  The  trouble 
with  this  method  is  that  it  too  only  works  for  linear  constraints.  Colors 
would  have  to  be  spaced  in  the  red,  green,  and  blue  coordinate  system,  which 
has  linear  boundaries,  and  then  transformed  into  L*u*v*  coordinates  to  cal¬ 
culate  the  distances,  much  in  the  manner  that  Carter  and  Carter's  algorithm 
does. 

Another  interesting  improvement  would  be  to  use  a  concept  called  "simulated 
annealing"  to  help  space  the  points.  In  this  method,  the  points  would  be 
able  to  actually  move  closer  together  in  order  to  eventually  achieve  optimal 
spacing.  Refer  to  Kirkpatrick,  Gelatt,  and  Vecchi  (1983)  for  more  on  this 
possibility. 

Other  less  drastic  refinements  might  be  to  make  the  program  interactive, 
more  user-friendly,  or  perhaps  include  a  color-naming  subroutine  that  would 
give  the  user  descriptive  names  of  the  colors  selected.  Carter  and  Carter's 
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program  includes  this  last  feature.  Also,  improvements  might  be  made  to  the 
speed  and  efficiency  of  the  program.  Consider,  for  example,  a  three-phase 
system  of  varying  the  maximum  and  minimum  step-sizes  and  the  step-size 
reducing  scheme  to  possibly  obtain  better  answers  more  quickly.  It  would 
consist  of  starting  out  with  large  maximum  and  minimum  step-sizes  and 
halving  the  step-size  at  each  iteration.  This  would  place  the  points  in  a 
roughly  spaced  configuration.  The  next  phase  would  have  a  much  smaller 
maximum  step-size,  reducing  by  four-fifths  to  the  same  minimum  step-size. 
This  would  refine  the  points  to  what  one  would  hope  to  be  a  rough  global - 
optimum  location  scheme.  Phase  3  would  concentrate  on  further  refining  the 
locations  and  the  accuracy  with  a  very  small  maximum  step-size  and  a  minus¬ 
cule  minimum  step-size.  The  composite  result  might  be  impressive. 

Yet  another  improvement  would  be  to  selectively  place  the  points  in  the 
corners  of  the  region  and  then  place  any  remaining  points  randomly  on  the 
faces,  rather  than  randomly  placing  all  points  in  the  interior  of  the  region 
as  an  initial  location  scheme.  This  procedure  would  conform  more  closely  to 
the  observations  made  in  Section  IV  regarding  simple  problems. 

It  is  recommended  that  future  efforts  include  extensive  statistical  tests 
comparing  the  algorithms  and  exploring  their  sensitivities  to  parameters. 

The  topic  is  interesting,  important,  and  satisfying,  and  anyone  else  who 
pursues  research  in  the  area  will  surely  enjoy  it. 
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Appendix  A 

COLOR  TRANSFORMATIONS 


This  appendix  briefly  describes  the  transformations  from  red,  green,  and 
blue  color  coordinates  to  the  CIE  L*u*v*  color  space.  If  we  let  (Y[^,  Yg, 

Yg)  be  the  luminances  for  the  CRT's  red,  green,  and  blue  guns,  respectively, 
which  produce  a  given  color,  then  the  first  task  is  to  find  something  called 
the  color's  tristimulus  values  (Xj,  Yj,  Zy).  This  transformation  is  given 
by  the  following  matrix  operation: 


'Xt' 

’>^R/yR 

’‘G/yc  xfi/yB 

"yr" 

Yt 

= 

1 

1  1 

Yg 

h 

2R/yR 

^G^yG  ^B/ye 

Yb 

where 

(xg,  yg,  zg)  =  the  chromaticity  coordinates  of  the  red  gun 

(xq,  yg,  zg)  =  the  chromaticity  coordinates  of  the  green  gun 

(xg,  yg,  Zg)  =  the  chromaticity  coordinates  of  the  blue  gun 

These  chromaticity  coordinates  correspond  to  the  locations  on  the  1931  CIE 

chromaticity  diagram  of  the  red,  green,  and  blue  colors  corresponding  to  the 
cathode-ray  tube's  guns,  and  are  values  between  zero  and  one.  In  addition, 

X  +  y  +  z  =  1  for  each  gun. 

The  next  step  is  to  find  the  u'  and  v'  values  for  the  color.  These  values 
correspond  to  the  color's  location  on  the  CIE-UCS  chromaticity  diagram  shown 
in  Figure  2.  The  transformations  are  given  by: 

u"  =  4x/(-2x  +  12y  +  3) 

V'  =  9y/(-2x  +  12y  +  3) 


where 
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X  =  Xj/(X-[-  +  Yj  +  Zy) 
y  =  Yy/(Xy  +  Yy  +  Zy) 


The  final  step  is  to  convert  from  (Yy,  u',  v')  to  (L*,  u*.  v*)  coor¬ 
dinates.  This  is  accomplished  by  the  following  set  of  equations: 

L*  =  116  (Yy/Yo)^''^  -  16 

u*  =  13  L*  (u'  -  Up) 

V*  =  13  L*  (V'  -  v^) 

where 

Yq  =  the  value  of  Yy  when  all  three  guns  are  set  to  their  maximum 
luminances,  producing  a  pure  white  color 

Uq  =  the  value  of  u'  for  a  white  corresponding  to  CIE  standard 
illuminant  Dg5  (=  .1954) 

vX  =  the  value  of  for  a  0,^  white  (=  .4697) 

U  DD 

The  definitions  for  Yg,  Ug,  and  v'  are  not  universally  agreed  upon  for  cases 
involving  self-luminous  displays.  The  convention  used  here  is  the  one  which 
Carter  and  Carter  used  and  has  been  chosen  because  it  enables  comparison  of 
the  algorithm's  results.  The  reverse  transformations  are  contained  in 
Appendix  B  (which  is  a  listing  of  the  program)  among  the  coding  of  sub¬ 
routine  "RG80UT." 
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Appendix 

FORTRAN  CODE  AND  DOCUMENTATION 


**t1liM*Mti*t*.*t***t**t**yj‘if '^***  *******  if  ill.  t**  *  t 

c 

i:  THF!;[S  FKOGRArt 

C  BY 

C  BOSS  KOLEY 

G0R-84B 


[tttt**t*Mt■t*.'l^*%t%t**1^*t*■^^*t.t*t.t*t***********t**t^^*1^.•»^**^^^^*t*■^*1^^* 

MAIN  PROGRAM 

THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  MAXIMIZE  THE  MINIMUM 
DIGTANCE  AMONG  N  POINTS  IN  A  CONVEX  THREE-DIMENSIONAL  REGION 
AND  THEN  SUCCESSIVELY  MAXIMIZE  THE  SECOND  SMALLEST  DISTANCE, 

THE  THIRD  SMALLEST,  AND  SO  ON,  UNTIL  ALL  POINTS  ARE  OPTIMALLY 
SPACED.  THE  REQUIRED  INPUTS  FOR  THIS  PROGRAM  ARE  THE  NUMBER 
OF  POINTS  BEING  SPACED  <N),  THE  DESIRED  INITIAL  STEP-SIZE  (SIZE), 
THE  NUMBER  OF  CONSTRAINTS  (FACES)  ON  THE  REGION  (F),  AND  THE 
EQUATIONS  OF  THOSE  CONSTRAINTS  IN  THE  FORM  AX  +  BY  ♦  CZ  <  K  . 

THE  PROGRAM  THEN  OUTPUTS  THE  OPTIMAL  VALUES  OF  X,  Y,  AND  Z, 

AND  THE  FIRST  TEN  (OR  SOME  OTHER  NUMBER  CHOSEN  BY  THE  USER) 
MINIMUM  DISTANCES. 

NOTE?  THROUGHOUT  THIS  PROGRAM,  X,  Y,  AND  Z  ARE  USED  TO  DENOTE 
THE  LOCATIONS  OF  THE  POINTS  BEING  SPACED  AND  CORRESPOND  TO 
L»,  U*,  AND  V*  FOR  THE  COLOR  SPACING  PROBLEM,  THIS  SHOULD  NOT 
BE  CONFUSED  WITH  THE  TRISTIMULUS  VALUES  OF  A  COLOR  WHICH  ARE 
COMMONLY  REFERRED  TO  AS  CAPITAL  X,  Y,  AND  Z  IN  THE  TRADE.  IN 
THIS  PROGRAM,  THE  TRISTIMULUS  VALUES  WILL  BE  REFERRED  TO  AS 
XTOTAL,  YTOTAL,  AND  ZTOTAL. 


THERE  ARE  SEVERAL  NOTES  TO  THE  USER  IN  THF  PROGRAM  WHICH 
ARE  USED  TO  SIGNIFY  PLACES  WHERE  YOU  MAY  BE  REQUIRED  TO  MAKE 
CHANGES.  THEY  ARE  SUMMARIZED  HERE  FOR  YOUR  CONVENIENCE. 

IN  THE  MAIN  PROGRAM,  THE  USER  IS  RFQUIRLD  TO  SUPPLY  THF 
NUMBER  OF  POINTS  TO  BE  SPACFD,  THE  SEED,  THE  MAXIMUM  STEP -SIZE, 
THE  MINIMUM  STEP-SIZE,  AND  THE  STEP-SIZE  REDUCING  SCHEME.  IN 
SUBROUTINE  PLACE-POINTS,  HE  MUST  SUPPLY  THE  MAXIMUM  AND  MIN¬ 
IMUM  VALUES  ALLOWABLE  FOR  EACH  OF  THE  THREE  COORDINATES.  IN 
SUBROUTINE  FIND-DISTANCE,  THE  USER  MUST  CHOOSE  THE  NUMBER  OF 
DISTANCES  HE  WISHES  TO  APPEAR  ON  THE  PF(OGF;AM'S  OUTPUT, 

IN  SUBROUTINE  CHECK-POINT,  THE  USER  IS  REQUIRED  TO  SUPPLY 
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C  )i;:  NIJMHI'.K  Ur  t  ACI  S  lllr>r  VtUUNn  IHL  Kf-UIDN  ANII  tHf  lOUATfOM:, 

C  FOK  I  HOBF  FALLS  IN  THI:  FOliM  AX  F  bY  F  127  K.  I  UK  IMF 

C  COLOR  SPACING  PROBLF_M,  THE  NUMBER  OF  FACFS  SHOUL  B  USUAL  LY 

C  BE  SFUEN.  IN  ABBITION,  All  OTHER  VARIABLES  STAY  THF  SAME 

C  EXCEPT  FOR  THE  VALUES  OF  K;  K(l)  IS  THE  MAXIMUM  ALLOWABLE 

C  LUMINANCE  YOU  WISH  TO  SET  FOR  THE  REB  GUN,  K(?)  =  THE  MAXIMUM 

C  GREEN  GUN  LUMINANCE,  K(3)  =  THE  MAXIMUM  BLUE  GUN  LUMIUANCL, 

C  ANB  K(7)  the  MINIMUM  TOTAL  LUMINANCE  YOU  WISH  TO  PERMIT  (IT 

C  SHOULB  PE  ENTIRLB  AS  A  NEGATIVE  VALUE). 

C 

C  IN  SUBROIIMNE  FIX  rOINT  THE  ONLY  REOUIREB  CHANGE  IS  I  HE 

C  VALUE  OF  THE  MAXIMUM  STEP-SIZE. 

C 

C  FOK  COLOR  SPACING  PROBLEMS,  SUBROUTINE  RGBOOT  IS  VERY 

C  IMPORTANT.  THERE,  THE  USER  MUST  SUPPLY  THE  CHROMATICITY 

C  COORBINATES  OF  THE  GUNS,  THE  MAXIMIM  LUMINANCE  ALLOWABLE 

C  <L-NOT),  ANB  THE  VALUES  OF  U  PRIME  NOT  AND  V~PRIME  NOT.  THAT 

C  IS  THE  EXTENT  OF  THE  CHANGES  THAT  NEED  TO  BE  MADE  FOR  ANY 

0  NORMAL  USE  OF  THIS  PROGRAM. 

C 

C 

C1^t*********t**•li^t*****t1^**t**1^*ttM*t************************■**t1^***tt* 

C 

C  VARIABLE  DEFINITIONS 

C 

c  REAL  variables: 

C  X<I)=  A  ONE-DIMENSIONAL  ARRAY  CONTAINING  THE  X-COORDINATE  OF 

C  POINT  I 

C  Y<I)=  THE  Y-COORDINATE  OF  POINT  I 

C  Z<I)*  THE  Z-COORDINATE  OF  POINT  I 

C  D<I,J)=  THF  EUCLIDEAN  DISTANCE  BETWEEN  POINT  I  AND  POINT  J 
C  CAPD=  THE  CURRENT  MINIMUM  DISTANCE  BETWEEN  ALL  PAIRS  OF  POINTS 
C  NEUD=  THE  NFW  MINIMUM  DISTANCE 

C  SIZE=  THE  STEP-SIZE 

C  ENDXd)*  THE  X-COORDINATE  OF  ENDPOINT  I,  WHERE  THE  ENDPOINTS 
C  ARE  THE  TWO  CLOSEST  POINTS  IN  THE  REGION 

C  ENDY<I)=  THE  Y-COORDINATE  OF  ENDPOINT  I 

C  ENDZ(I)=  THE  Z-COORDINATE  OF  ENDPOINT  I 

C  AL.TX<I,J)=  THE  JTH  ALTERNATE  LOCATION  FOR  THE  X-COORDINATt  OF 

C  ENDPOINT  I 

C  ALTY<I,J)=  THE  JTH  ALTERNATE  LOCATION  FOR  THE  Y-COORDINATE  OF 

C  ENDPOINT  I 

C  ALTZ<I,J)=  THF  JTH  ALTERNATE  LOCATION  FOR  THE  Z-COORDINATE  OF 

C  ENDPOINT  I 

C  MEMX(I)=^  THE  X-COORDINATE  OF  POINT  I  FROM  THE  PREVIOUS  ITERATION 

C  MEMY(I)=  THE  Y-COORDINATE  OF  POINT  I  FROM  THE  PREVIOUS  ITERATION 

C  MEMZ(I)=  THE  Z-COORDINATE  OF  POINT  I  FROM  THE  PREVIOUS  ITERATION 

C  DUM1,DUM2,DUM3=  DUMMY  VARIABLES  REQUIRED  FOR  SUBROUTINE  RGBOUT 

C  BUT  NOT  USED  IN  THE  MAIN  PROGRAM 

C  SFED=  THE  STARTING  VALUE  FOR  THE  RANDOM  NUMBER  GENERATING  SEQUENCE 
C  IN  THE  CYBER  KNOWN  AS  SUBROUTINE  'RANSET' 

C 
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nooi-  nnnnnnnnnn 


i:  fNTt  Of:K  '.Mlslr’^hLLij  : 

i:  i:it;  mumblr  or-  points  PLOuiurNC  sp^cinc: 

C  PNDFTl^  THE  NUMBER  ASSOCIrtlED  WIIH  ENUPOJNT  1 

r:  ENIH'TI’^  THE  NUMBER  ASSOCIATED  UITH  EINDPOrNT  I’’ 

C 

C  CHARACTTR  VARIABLES: 

C  STATUS--  THE  STATUS  OE  THE  PROGRAM;  'HUN'  IF  THE  STOPPINH 

C'  CR  HER  ION  HAS  NOT  BEEN  MET,  AND  'DONE'  IF  IT  HAS 

L  STATPT(I):^  IHir  SIATUS  OF  POINT  i;  'F'  FOR  FIXED,  AND  'U' 

C  FUR  UNI  I, VP D 

C 

C  EXTERNAL  VARIABLES: 

C  RANI-:-:  VARIABLE  USED  IN  THE  RANDOM  NUMBER  OENERATING  SUBROUTINE 

C  IN  I  HE  CYBER  CALLED  'RANSET' 

C 

C 

Ct*tyif***t*****t*t'i(***************t*tt**.i(*yift****t**t********tt*.tt**t**t 

c 

C 

PROGRAM  SPACPT 

REAL  X(50) ,Y(50) ,Z<50),D<50,S0),CAPD,NEWD,SI2E,ENnX(2) ,ENDY(2>  , 

I  ENDZ(2) ,ALTX( 2,27) ,ALTY( 2,27) ,ALTZ( 2,27) ,MEMX (50) ,MEMY(50) , 

S  Mr MZ(50) ,DUM1<50) ,DUM2<50) ,DUM3(50) ,SEED 

INTEGER  N,ENriPTl,FNDPT2 
CHARACTER  STATUS».4 ,  STATPT  ( 50  )  *1 
EXTERNAL  RANF 
OPEN  <3,FILE*'SPC0UT' ) 

REWIND  3 
10*3 
100*3 

—  IMPORTAHT  NOTE  TO  USER  — 

THE  USER  MUST  SUPPLY  THE  VALUE  OF  N  AND  SIZE  HERE,  AND  THE  SIZE 
AGAIN  IN  SUBROUTINE  'FIX-POINT. '  I  RECOMMEND  USING  A  STEP-SIZF 
APPROXIMATELY  EQUAL  TO  ONE-HALF  THE  DISTANCE  BETWEEN  THE  TWO 
FARTHEST  POINTS  IN  THE  REGION,  BUT  THE  USER  IS  ENCOURAGED  TO  TRY 
THE  SAME  problem  WITH  SEVERAL  DIFFERENT  MAXIMUM  STEP-SIZES  TO 
INCREASE  THE  CHANCES  OF  GETTING  AN  OPTIMAL  SOLUTION. 

SIZE*130.0 
N=10 

STATUS  -- '  DONE  ' 

DO  100  |:vl,N 

STATPT ( I )='U' 

MEMX(I)*0.0 
MEMY(  I):^0.0 
MEMZd  )*0.0 
CONTINUE 

—  IMPORTANT  NOTE  TO  USER  — 
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n  n  n  o 


c  !i;i  M  "  ‘,'t^l  iH.  Kll-  llll  Ini'lnih  .liuli't.  I  :i  NI-,  I  r  f  ( !( '. 

L  scriiii:  ni:e  wh/i.ii  i ..'imliuhly  CHOiists  i;il  iwiTirti  locaiilins  Mjiv  nu,. 

c;  l-■(J[Nls.  ir  rr.  i-:i- roMMi-Nruiii  rur-.i  h'ji  iii  it  ftUNb  kh  pi  ismihmi  n  roi; 

c  1  f'IsOf'Lt  MS,  rN  UHTCH  Crt'bP:  IHt  Un^H  SHQUl  Ii  si  1  D  FUK 

C  I  YiCI-l  ruiN.  ir  I'r’iP'  KF  rtNY  NfNK  D  I  (i  I  T  iNriljlR. 

L' 

SI  P  1.1-^  7(r.''.>is  1  :*i  ,0 
wKJTFCi,  1  iO)sr::D 

1 10  p  ORMA  r  ( 1  sv ,  ■ '  I  II  sf- 1- 11  •  ,p  1 1 ,  n 

i:  Mi/s  srrcTf'11,1  c.'ii  i  t"  ri ic s  rur  /Mriioi.  locop  cons  ’  ok  mf  fopnps  ond 

C  FORHAr:  ni''Ti  PR  IN  IS  1 1  IT  I  GCA  T  1  CiNS  . 

CALL  RANSCr  1  SI;  FD) 

CAI  L  l-'l  AC  P  T  (  N  ,  ,  Y  v  /  ,  SPl.H  H  ) 

UR  [  I  f.  (  S  ,  1  00  ) 

]  00  FORMAT  <  n  '  ITILSE  hRF  THE  STAR  I  INI.  I.UCAT  IONS  '  ) 

CALL.  RR[NT(X,Y,7,  ,N) 

CALL  RGROU  1  (  X  v  Y  ,  7  ,  N  ,  HI  INI  ,  DUMO  ,  DlJM.i  ,S  r  A  I  US  ) 

STATUS^ 'RUN' 

C  THE  MEAT  OF  THE  l-RGGRAM  BEGINS  HERE 

130  CALL  F IMHD  (  X  ,  Y  ,  Z  ,  N  ,  STATf'T  ,  CARD ,  ENDX  ,  ENDY  ,  ENDZ  ,  ENDPT 1 ,  ENDPT2 , 

X  STATUS) 

CALL  DEI- ALT  (ENDX,  ENDY,  ENDZ,  SIZE,  ENDPT  1  ,ENDPT2,STATPT  , 

8  ALTX,ALTY,ALr2) 

CALL  CHKALT  (  X  ,  Y  ,  Z  ,  N  -  •'APD  ,  ALTX  ,  ALTY  ,  ALT7  ,  STATPT  ,  ENDPT  1  ,  ENDPT2  ,  NEUD  ) 

C  IF  THE  MINIMUM  DISTANCE  IS  GREATER,  THEN  KEEP  ITERATING 
IF  <NEWD  .GT.  CAPD)  THEN 
CAPD=NEWD 
GO  TO  130 

C  IF  NOT,  THEN  REDUCE  THE  STEP-SIZE 
ELSE 

SIZE-SI2I-/2.0 
IF  (SIZE  .LE.  .0001)  THEN 

NOTE?  THE  USER  MAY  WISH  TO  USE  A  DIFFERENT  MINIMUM  STEP-SIZE. 

FOR  THREE  DECIMAL  PLACES  OF  ACCURACY,  T  USE  .0001.  YOU 
MAY  ALSO  WISH  TO  USE  A  SCHEME  OTHER  THAN  HAl.UING  THE  STEP- 
SIZE,  LIKE  REDUCING  IT  BY  THF^EE-FOURTHS  INSTEAD. 

CALL  FIXPT<X,Y,Z.N,STATPT,ENDPT1 ,ENDPT2, SIZE, STATUS. 

8  MEMX,MEMY,MEMZ) 

IF  (STATUS  .EQ.  'RUN')  THEN 
GO  TO  130 

ELSE  IF  (STATUS  .EO.  'DONE')  THEN 
GO  TO  140 
END  IF 

ELSE 

GO  TO  130 
END  IF 
END  IF 

C  THIS  SECTION  CONTROLS  THE  FORMAT  OE  THE,  RESULTS 
140  URITE(3,150) 

150  FORMAT (15X, 'THIS  IS  THE  FINAL  ANSWER') 

CALL  PRINT(X,Y,Z,N) 

CALL  RGBOUT ( X , Y , Z , N , DUMl , PUM2 , DUM3 , STATUS ) 
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Lrti.  I  I  lNIiH  <  ^  ,  f  .  c  ,  H  >  rp  I  ,rnK-\i  ,r  Nd  <  ,  i  I'/OY  ,  ;  rUi/  ,  i  I'lfir  I  i  ,  ;  i  M'.  '  !  , 

X  sihri's) 
tND 
C 
C 

C 

C  V'.HEiROUTlNf-  PI..ACF.  PO(N!;> 

C 

c  Till  i'Iuttosi::  ni-  iniy  haby  is  to  RnNOom  y  iT.rtOfc'  thiv  PH  imp; 

C  IN  THE  REGION.  THE  LITTLE  TYKE  fVECEIL'ES  THE  NIJMHER  OF  POINTS 

C  (Nl  AS  AN  INITJT,  ANH  OORGLES  UP  TUP  STARTING  l  OGATIDNS  TOR  Till 

C  POINTS  (X,  Y,  AMU  Z)  AS  OUTPUTS. 

C 

c  note:  this  subroutine  uses  the  random  number  generating  function 

C  IN  THE  CYBER  CALLED  'RANSET.'  MOST  MAIN-PRAME  COMPUTERS  HAVE 

C  A  RANDOM  NUMBER  GENERATING  CAPABILITY,  SO  CONSULT  YOUR  LOCAL 

C  COMPUTER  EXPERT  FOR  DIRECTIONS  ON  HOW  TO  USE  IT  ON  YOUR  FAVORITE 

C  SYSTEM. 

C 

C 

C*tt**t****'***********tt*t*****t***t*********************************** 

C 

c  VARIABLE  DEFINITIONS 

C 

c  REAL  variables: 

C  X<I)=  THE  INITIAL  X-COORDINATE  OF  POINT  I 

C  Y(I)=  THE  INITIAL  Y-COORDINATE  OF  POINT  I 

C  Z<I)=  THE  INITIAL  Z-COORDINATE  OF  POINT  I 

C  YRED(I)=  THE  RED  PHOSPHOR  LUMINANCE  FOR  POINT  I 


C  YGREEN<I)=  THE  GREEN  PHOSPHOR  LUMINANCE  FDR  POINT  I 

C  YBLUE(I)=  THE  BLUE  PHOSPHOR  LUMINANCE  FOR  POINT  I 

C  SEED=  THE  STARTING  VALUE  FOR  THE  RANDOM  NUMBER  GENERATING  SEOUENCE 

C  IN  THE  CYBER  KNOWN  AS  SUBROUTINE  'RANSET' 

C 

C  INTEGER  variables; 

C  N=  THE  NUMBER  OF  POINTS  BEING  LOCATED 

C  J=  TI-IF  NUMBER  OF  POINTS  BEING  TRANSFORMED  BY  SUBROUTINF  RGBOUT' 

C 

C  CHARACTER  VARIABLES: 

C  POINT=  THE  LOCATION  OF  THE  POINT;  EITHER  'IN'  OR  'OUT'  OF  THE 

C  REGION 

C  DUM1=  A  DUMMY  VARIABLE  REOUIRED  FOR  SUBROUTINF  'RGBOUT'  BUT  NOT 

C  USED  IN  THIS  SUBROUTINE 

C 

C  EXTERNAL  VARIABLES: 

C  RANF=  A  VARIABLE  USED  IN  CONJUNCTION  WITH  'RANSET'  TU  GENERATE 

C  A  RANDOM  NUMBER  BETWEEN  ZERO  AND  ONE 

C 
C 

C 
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SUBl>niJ  T  J  Mb!  b  l.ACPT  (  N  ,  X  ,  i  ,  Z  ,  ) 

REAL  X(50)  ,  Y(50)  ,/(50)  ,SEED,YRt:ii(50)  ,  YGFtE  EN  <  ^,0  )  ,  YruUf  (  bO  ) 

INTEGER  N,J 

CHARACTER  POINT *1 , DUM1»4 
EXTERNAL.  RANK 
no  TAG  I=1,N 

IZO  X(  I )  =:RANF(  )«<  100.0-9.0)  f9.0 

Y  (  I  )=I<ANF  <  )*(  175.0  +  75.0  )-7fj.0 
7( I )^RANF< )»( 100. 0+1 25.0 >-125.0 

--  IMPORTAMV  NOTE  TO  USER  -- 

THE  VARIARLE  'RANE'  GIVES  VALUES  IN  THE  RANGE  FROM  ZERO  TO  (TNIi  . 

IF  THE  USER  WANTS  VALUES  OUTSIDE  OF  THAT  RANGE,  THE  EQUATIONS 
FOR  X,  Y,  AND  Z,  SHOULD  LOOK  LIKE  THIS: 

X<I)=RANF< )*<XHAX-XMIN)+XMIN 
Y<  I  )=RANF(  )».<YMAX-YMIN)+YMIN 
Z<I)=RANF( )»(ZMAX-2MIN)+ZMIN 

WHERE  MAX  AND  MIN  ARE  THE  MAXIMUM  AND  MINIMUM  POSSIBLE  VALUES 
FOR  X,  Y,  AND  2,  RESPECTIVELY. 

J=1 

CALL  F;GB0UT(X<1),Y<I>,Z<I),J,YRED(I),YGREEN(I),YBLUE(I) ,DUM1) 
CALL  CHEKPT ( YRED< I ) , YGREEN< I ) , YPLUE< I ) , POINT) 

IF  (POINT  .EQ.  'OUT')  THEN 
GO  TO  170 
END  IF 

60  CONTINUE 
END 


SUBROUTINE  FIND-DISTANCE 


THF.  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  FIND  THE  MINIMUM  DIS¬ 
TANCE  BETWEEN  AI.L  N(N-l)/2  PAIRS  OF  (UNFIXED)  POINTS.  THE 
INPUTS  FOR  THIS  PROGRAM  ARE  THE  LOCATIONS  OF  THE  POINTS  (X, 

Y,  AND  Z),  THE  STATUS  OF  THOSE  POINTS  AS  TO  WHETHER  THEY  ARE 
FIXED  OR  UNFIXED  (STATPT),  AND  THE  NUMBER  OF  POINTS  (N).  THE 
PROGRAM  THEN  OUTPUTS  THE  MINIMUM  EUCLIDEAN  DISTANCE  BETWEEN 
ALL  PAIRS  OF  (UNFIXED)  POINTS,  THE  TWO  POINTS  THAT  CORRESPOND 
TO  THAT  MINIMUM  DISTANCE  (ENDPTl  AND  EN0PT2),  AND  THE 
COORDINATES  OF  THOSE  POINTS  <ENDX,  ENDY,  AND  ENDZ). 


VARIABLE  DEFINITIONS 


REAL  variables: 
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i:  X'.''-  iMf:  X  nr  iiiinr  r. 

c  Y<[)^  nil-;  Y  cnnnriiNATR  or  i-oini  i 

.c  z<n^  the:  z-cooHniNAiF.  of  point  i 

C  CAPH-  Tlir.  MINIMUM  EUri  IHFAN  DtyTMNCE;  HHIUEEN  MIL  (UNFIXFIO  PAIRS 

C  OF  POINTS 

c  i:NriX(i)=  THi-  y -cnriRniNATL-  or-  entipoint  i 

r:  ENUYd)  ^  T  HE  Y-COORU  [NATE  OF  ENDPOINT  I 

C  FNn/dT^^  THE  /-COORDINATE  OP  ENDPOINT  I 

C  D([,J)=  THE  euclidean  DISTANCE  FROM  POINT  I  TO  POINT  J 

C  ITISd'n-"  IHE  MU  SMALLEST  UlSTANCI; 

C 

C  INTEGER  UARIAHLFSr 

C  N=^  THE  NIJMPER  OF  POINTS 

C  i:NliPTl=^  THF  NUMDFR  ASSOCIATED  WITH  ENDPOINT  1 

C  ENliPT2^  THE  NIJMDER  ASSOCIATED  WITH  ENDPOINT  2 

C  K=  THF  COUNTER  USED  IN  ORDERING  THE  DISTANCES 

C  L.  =  THE  NUMBER  OF  UISFANCES  THE  USER  WISHES  TO  PRINT 

C  EYE<K)=  AN  ENDPOINT  COORESPONDING  TO  THE  KTH  SMALLEST  Dd,J) 

C  JAY(K)-^  THE  OTHER  ENDPOINT  COORESPONDING  TO  THE  KTH  SMALLEST  Dd,J) 

C 

C  CHARACTER  VARIABLES: 

C  STATPT<I)=  THE  STATUS  OF  POINT  i;  EITHER  FIXED  OR  UNFIXED 

C  STATUS=  THE  STATUS  OF  THE  PROGRAM;  EITHER  'DONE'  OR  'RUN' 

C 

C 

Ct$*****t**f***************************$*************$**************** 

c 

c 

SUBROUTINE  FINDD < X , Y , 7 , N , STATPT , CAPD r ENDX , ENDY , ENDZ , ENDPT 1 , ENDPT2 , 
t  STATUS) 

REAL  X<S<>) ,Y<50) ,Z<50) ,CAPD,ENDX(2) fENDY<2) fENDZ<2) ,0(50,50) , 

X  DIS<300) 

INTEGER  N , ENDPT 1 , ENDPT2 , K , L , EYE ( 300 ) , JAY ( 300 ) 
character  STATPT  <  50 )  *1 ,  STATUSM 
CAPD= 10000.0 
DO  200  1=1, N-1 

DO  210  J=:i  +  ],N 

C  UNLESS  BOTH  POINTS  ARE  FIXED,  CALCULATE  THE  DISTANCE  BETWEEN  THFM 
IF  ((STATPTd)  .EQ.  'U')  .OR.  (STATPTCJ)  .EQ.  'U')) 

S  THEN 

C  THIS  IS  THE  EUCLIDEAN  DISTANCE  FORMULA 

D  d , J ) =SQRT ( ( X  d ) -X ( J ) ) »»2+ ( Y ( I ) -Y ( J ) ) »*2+ 

X  <Zd)-Z<  J)  X-.*!') 

IF  <Dd,J)  .LE.  CAPD)  THEN 
CAPD=Dd,J) 

FNDPTl=t 
ENDPT2=J 
END  IF 
END  IF 

210  CONTINUE 

200  CONTINUE 

ENDX(l)-^X<ENIiPTl) 
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FNIIY(  1  )  -  r  (|;NMPI  !  ) 
t:Mr.i2(  1. )  2  a  Ntipi  i ) 

ENIiX(21--X(I!:NDIT?) 

r:MtiY(;:)^Y(n;MiipT2) 

F.N1IZ(:>)-7(PNDPT2) 

c  THIS  si:t:iiON  npriKus  the:  uisrANCEs  and  pi;inis  the  l  smallest 

L  DISTANCES  TO  OUTPUT  WHEN  THE  PROGF'iAM  IS  DONE 
;[E  (STATUS  .i:0.  DONE")  THEN 
L  =  10 

C  THE  USER  CONTROLS  THE  NUMBER  Of-  DISTANCES  HE  W:i.SHE5  TO  Sl-E  BY 
C  CHANGING  THE  VALUE.  OF  L  HERE,  HE  CAN  CHOOSE  TO  SEE  UP  TO  ;300 
C  DISTANCES. 

DO  220  K- 1,300 

£irS(K) -1000.0 
220  CONTINUE 

DO  230  1=1, N-1 

DO  240  .1=1  +  1, N 

DO  250  K--l,l.. 

IF  <D(I,J)  .LT.  DIS(K))  THEN 
DO  260  M=L,K+1,-1 

DIS(M)=DIS(M-1) 

EYE(M>=EYE(M-1> 

JAY(M)=.JAY(M-1) 

260  CONTINUE 

DIS<K>=Ii(I,  J) 

EYE(K)=I 
JAY(K)=J 
GO  TO  240 
END  IF 

250  CONTINUE 

240  CONTINUE 

230  CONTINUE 

WRITE(3,280)L 
DO  290  K=1,L 

URI TE ( 3 , 295 ) EYE  <  K ) , JA Y ( K ) , D I S ( K ) 

290  CONTINUE 

URITE<3,270) 

270  FORMAT ('  ') 

280  F0RMAT(9X, 'THESE  ARE  THE  ',13,'  SMALLEST  DISTANCES') 

295  FORMAT! 19X, 'D< ' , 12 , ' , ' , 12 , ' ) = ' , F7 . 3 ) 

END  IF 
END 
C 
C 

C 

C  SUBROUTINE  DEFINE-ALTERNATIVES 

C 

C  THE  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  DEFINE  THE  VARIOUS 

C  ALTERNATIVE  LOCATIONS  FOR  ENDPOINTS  1  AND  2.  THE  PROGRAM 
C  IS  PASSED  THF  COORDINATES  OF  THE  ENDPOINTS  AND  THE  STEP-SIZE 

C  AND  OUTPUTS  A  TWO  BY  TWENTY-SEVEN  ARRAY  CONTAINING  THE 
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CUORII  (NrtTt'j  or  IHK  Al  ,TI- h:NA  n  VE  S  .  Illl-rih:  Al.  1 1.  RN(M  [(JL  I  II 

(JN  A  fUff.  WHurjF.  aiDfc:s  are  twice  thl  :;if.p-size  in  niciANCE 

WITH  THE  ENIIPOINT  AT  IT!J  CENTER. 


C 
C 
C 
C 
C 

Ct********t*'*^*'‘r-1f*'V*t**MtM*ft*t***1l-****tt*-1tittt*ttt**t****t**1li1ti**ttt*t** 

c 

C  'VARIABLE  DEEINITIUNy 

C 

C  REAL  yARIABI  Eij-r 

C  ENDX<I)=  THE  X -COORDINATE  OF  ENDPOINT  I 

C  ENDY(I)-  THE  Y  COORDINATE  OF  ENDPOINT  I 

C  END7(I)=  THE  Z-COORDTNATE  OF  ENDPOINT  I 

C  SIZE=  THE  STEP-SIZE 


c 

ALTXd 

,J)=  THE 

X -COORDINATE 

OF 

THE 

JTH 

ALTERNATIVE 

FOR 

ENDPOINT 

I 

c 

ALTYd 

,J)=  THE 

Y-COORDINATE 

OF 

THE 

JTH 

ALTERNATIVE 

FOR 

ENDPOINT 

I 

c 

ALTZd 

,.!)=  THE 

Z-COORDINATE 

OF 

THE 

JTH 

ALTERNATIVE 

FOR 

ENDPOINT 

I 

C 

C  INTEGER  variables: 

C  ENDPT1=  THE  NUMBER  CORRESPONDING  TO  ENDPOINT  1 

C  ENDPT2*  THE  NUMBER  CORRESPONDING  TO  ENDPOINT  2 

C  P*  THE  INDEX  FOR  THE  ALTERNATE  LOCATIONS 

C 

C  CHARACTER  VARIABLES: 

C  STATPT<I)=  THE  STATUS  OF  POINT  i;  EITHER  FIXED  OR  UNFIXED 

C 

C 

»»«««««««»««««««»«»«««««««« 

c 

c 

SUBROUTINE  DEFALT ( ENDX , ENDY ,END2 , SIZE f ENDPT 1 , ENDPT2 , STATPT , ALTX , 

S  ALTY,ALTZ) 

REAL  ENDX (2) ,ENDY<2) ,ENDZ<2) ,SIZE,ALTX(2,27) ,ALTY(2,27) , 

X  ALTZ<2,27) 

INTEGER  ENDPTl,r''DPT2rP 
CHARACTER  STATPl , 50 ) *1 

C  THESE  LOOPS  LOAD  THE  ALTERNATE  LOCATIONS  FOR  ENDPOINTS  1  AND  2 
DO  300  T=l,2 
P=0 

DD  310  S=-l,l 

DO  320  R*-l,l 
DO  330 

P=P+1 

ALTX(T,P)=ENDX(T)-S*SIZE 
ALTY<T,P)=ENDY<T)-R»SIZE 
ALTZ<T,P)=ENDZ<T)-Q*SIZE 
330  CONTINUE 

320  CONTINUE 

310  CONTINUE 

300  CONTINUE 

C  IF  ONE  OF  THE  ENDPOINTS  IS  FIXFD,  THEN  IT  HAS  NO  ALTERNATE  POSITIONS 
IF  <STATPT(EN0PT1)  .EQ.  'F')  THEN 


n  n  n  n 


i.m  5^0 

ALTX(  l  ,S)  1:;NDX(  1  ) 

ALTY<l,S)-ENIiY(l) 

ALrZ(  1  ,S)-ENtiZ(  1) 

340  CONTINUE 

ENU  II-' 

IF  (statpt<fnuf-t::')  .eq.  'F')  then 

no  3“30  S  =  l,27 

At  TX(2,S>=ENr.iX<2' 

AL  rY(2,rO-FNtiY(2) 

Al.TZ(2,S)==FNliZ(2) 

350  CONTINUE 

END  II- 
END 
C 
C 

Cilliilit**tt1lit**tt*t'llit*t******'ttM*tt*tt*'lllittt***tt*t*t*****t*ilit*t***t**ttt*t 

C 

C  SUBROUTINE  CHECK-ALTERNATIVES 

C 

C  THIS  SUBROUTINE  CHECKS  THE  VARIOUS  ALTERNATIVES  TO  SEE 

C  IF  THEY  LIE  IN  THE  REGION,  THEN  TESTS  THE  REMAINING  CANDIDATES 

C  TO  SF.E  IF  THEY  INCREASE  OR  DECREASE  THE  DISTANCE  BETWEEN  THE 

C  ENDPOINTS.  IF  THE  DISTANCE  IS  INCREASED,  THIS  PROGRAM  CALLS 
C  SUBROUTINE  'FIND-DISTANCE'  TO  SEE  IF  THE  OVERALL  MINIMUM 
C  DISTANCE  IS  INCREASED  OR  DECREASED.  IF  SO,  THE  ENDPOINTS 
C  ARE  MOVED  TO  THE  ALTERNATE  LOCATION  THAT  MOST  INCREASES  THE 
C  OVERALL  MINIMUM  DISTANCE  WITH  TIES  GOING  TO  THE  ALTERNATIVE 

C  THAT  INCREASES  THE  ENDPOINT  DISTANCE  THE  MOST.  IT  THEN  PASSES 

C  THOSE  NEW  LOCATIONS  TO  THE  MAIN  PROGRAM. 

C 

C 

C*t**tt*****t******1ti****************M***************f-***************** 

VARIABLE  DEFINITIONS 

REAL  variables; 

C  X(I)=  THE  X-COOROINATE  OF  POINT  I 

C  Y<I)=  THE  Y-COORDINATE  OF  POINT  I 

C  Z<I)-  THE  Z-COOROINATE  OF  POINT  I 

C  CAPD=^  THE  MINIMUM  DISTANCE  BETWEEN  ALL  PAIRS  OF  (UNFIXED)  POINTS 

C  ALTX<I,J)=  THE  X-COOROINATE  OF  THE  JTH  ALTERNATIVE  FOR  ENDPOINT  I 

C  ALTY<I,J)=  THE  Y-COORDINATE  OF  THF  JTH  ALTERNATIVE  FOR  ENDPOINT  I 

C  AI..TZ<I,J)=  THE  Z-COOROINATE  OF  THE  JTH  ALTERNATIVE  FOR  ENDPOINT  I 

C  NEWD^  THE  NFW  MINIMUM  DISTANCE  BETWEEN  ALL  THE  POINTS 

C  CHKX<I)=  THE  X-COORDINATE  OF  POINT  I  BEING  CHECKED  AS  A  NEW 

C  LOCATION  SCHEME  BY  SUBROUTINE  'FIND-DISTANCE' 

C  CHKY(I)=  THE  Y-COORDINATE  OF  POINT  I  BEING  CHECKED 

C  CHKZ(I)=  THE  Z-COORDINATE  OF  POINT  I  BEING  CHECKED 

C  DUM1,DUM2,DUH3»  DUMMY  VARIABLES  THAT  ARE  RETURNED  BY  SUBROUTINE 
C  'FIND-DISTANCE'  BUT  ARE  NOT  USED  IN  THIS  SUBROUTINF 

C  MIND^  THE  MINIMUM  DISTANCE  CALCULATED  BY  SUBROUTINE  'FIND-DISTANCE' 


56 


c  HflH  THt:  LUCATION  SCHt  rtfc  DLI  TNI£lt  LO'  I'llKXrCHKi,  liNT  CHKZ 

C  MAXFND=  THE  VrtRIABlE  USED  rtS  A  nE-BKEAM^ft  FOF<  (=.L  lERfMTrv'ES  'JIIH 

C  IDENTICAL  MINIMUM  DISTANCES 

C  ENDPTD--=  THE  DISTANCE  BETWEEN  ENDPOINTS  1  AND  2 

C  YRED<I,.I)=  THE  VALUE  OF  ALTX(I,J)  IN  RE  D -GHFEN -BLUE  SPACE 

i:  YGRI-EN<  T  ,  J)-  THE  VAI.UE  OF  ALTYd.J)  IN  RGB  SPACE 

C  YDLIJl:  (I,J)  -  THE  VALUE  OF  AL!Z(I,J)  IN  RGB  SPACE 

C 

C  INTEGER  variables: 

C  ENDPT1--  THE  NUMBER  ASSOCIATED  WITH  ENDPOINT  1 

C  ENDPT2=  THE  NUMBER  ASSOCIATED  WITH  ENDPOINT  2 

C  N=  THE  NUMBER  OF  POINTS 

C  DUM4,DUM5^  NUT  USF.D  IN  THIS  SUBROUTINE 

C 

C  CHARACTER  VARIABLES: 

C  STATPT(I)=  THE  STATUS  OF  POINT  i;  EITHER  FIXED  OR  UNFIXED 

C  LOCALT(I,J)=  THE  LOCATION  OF  THE  JTH  ALTERNATIVE  FOR  ENDPOINT  i; 

C  EITHER  IN  OR  OUT 

C  DUM6=  A  VARIABLE  REQUIRED  FOR  SUBROUTINE  'RGBOUT'  BUT  NOT  USED 

C  IN  THIS  SUBROUTINE 

C 
C 

C 

C 

SUBROUTINE  CHKALT < X , Y ,Z ,N,CAPD,ALTX, ALTY ,ALT2 ,STATPT ,ENDPT1 , 

I  ENDPT2,NEWD) 

REAL  X(50),Y<S0),Z<50),CAPDfALTX<2,27),ALTY(2r27),ALTZ(2f27), 

S  NEWD,CHKX<50) ,CHKY<50) rCHKZISO) , DUMl ( 2) , DUM2( 2) ,0UM3<2> , 

S  MIND,HAXEND,ENLiPTDfYRED(2»27)  ,YGREEN(2,27)  rYBLUE<2f  27) 

INTEGER  ENDPTl ,ENDPT2,M,DUM4,DUM5 
CHARACTER  ST ATP  T  <  50 ) * 1 , LUCALT  <  2  r  27 ) »3 , DUM6*4 
NEUD=CAPO 
MAXEND=CAPD 
DO  400  1^1, N 

CHKX<I)=X<I) 

CHKY<I)=Y<I) 

CHKZ<I)=Z<I) 

400  CONTINUE 

C  THIS  LOOP  CHECKS  ALL  54  ALTERNATIVES  TO  SEE  IF  THEY  ARE  IN  OR  OUT 
C  OF  THE  REGION 

DO  410  S=l,2 

DO  420  T=l,27 

CALL  RGBOU  T  <  A!  TX  <  S  ,  T  )  ,  Al.TY  <  S ,  T  )  ,  ALTZ  ( 5 ,  T )  ,  1  r  YRED <  S ,  T  )  , 
S  YGREEN<S,T) ,YBLUE(S,T) ,DUM6) 

CALL  CHEKPT  <  YRED  <  S  r  T ) , YGREEN  <  S , T ) , YBLUE ( S , 1 ) » 

S  LOCALT(S,T)) 

420  CONTINUE 

410  CONTINUE 

DO  430  T=l,27 

IF  (L0CALT<1,T)  .EQ.  'IN')  THEN 
DO  440  S=l,27 
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:il-  (T  .Nh  .  S  .ANIi.  1.  DOM  I  ■:  :',S  )  .l  U.  '  IN  )  I HKN 
D  IK  fOTH  ENrif'O  [NTS  ARK  STILL  IN,  IHRN  I'ALCUI.  rtTL  IMG  NLW  DXSIAHL'L 
C  BETWEEN  THEM 

ENr»PTH-SQRT(  (ALIX(  i  ,  D  -Al  TX(2,r>) 

%  <  ALTY<  1 ,  T)  -ALTY(2,i:.)  )  **2+ (  ALTZ  (  1  ,  T  )  - 

5  ALTZ<2,'o)  )*»2) 

C  IF  THAT  DISTANCE  IS  GREATER,  THEN  CHECK  THE  OVERALL  MINJNUH  M I S I ANCE 

IF  (ENDPlIi  .GT.  NEUD)  THEN 
CIIKX(ENriPTl  )  rft|_  !X(U  T) 

CHKY(EMDPT1 )^ALTY(  1  ,  D 
CHKZCENDPTl  L^ALTZd  ,  I) 

CHKX  (  f:MIiPT2  )  =AL  I  X  (  2  ,  S ) 

CHKY  ( ENriPT2 )  =AI  T  Y  (  2  ,  S  ) 

CHKZ(ENIiPI2)-=AL  J  Z(2,S) 

CALL  FINDD(CHKX,CHKY  ,CHKZ,N,STATPTrMINri, 

6  DUM 1 ,  riUM2 ,  BUMS  ,  DUM  ■» ,  DUHS  ,  D1JM6  ) 

C  IF  THE  MINIMUM  DISTANCE  IS  THE  SAME  BUT  WITH  A  GREATER  ENDPOINT 
C  DISTANCE,  OR  IF  THE  MINIMUM  DISTANCE  IS  BETTER,  THEN  MOVE  THE 
C  ENDPOINTS 

IF  ((MIND  .GT.  NEWD)  .OR.  (MIND  .EO. 

X  NEWD  .AND.  ENDPTD  .GT.  MAXEND))  THFN 

NEUD-MIND 

MAXEND^ENDPTD 

X(ENDPT1>=ALTX(1,T) 

Y(ENDPT1)»ALTY(1,T) 

2(ENDPT1)-ALTZ(1,T) 

X(ENDPT2>=ALTX(2,S) 

Y(ENDPT2)=ALTY(2,S) 

Z(ENDPT2)=ALTZ(2,S) 

END  IF 
END  IF 
END  IF 

440  CONTINUE 

END  IF 

430  CONTINUE 
END 


*t-*t**t*t*****iii*tt*t*****t***M****ttt*tty^tt*^**tt**yii.tiii'iiii***tt.iiit.tM****t 

SUBROUTINE  CHECK-POINT 

THIS  SUBROUTINE  CHECKS  POINTS  TO  SEE  IF  THEY  ARE  IN  OR 
OUT  OF  THE  REGION.  IT  CHECKS  THEM  AGAINST  A  SET  OF  EQUATIONS 
(CONSTRAINTS)  OF  THE  FORM  AX  +  BY  +  CZ  ^  K  WHICH  MUST  BE 
SUPPLIED  BY  THE  USER,  AND  THEN  OUTPUTS  A  CHARACTER  VARIABLE 
WHICH  IDENTIFIES  WHETHER  THE  POINT  IS  'IN'  OR  'OUT.' 

NOTE:  THIS  SUBROUTINE  IS  FOR  LINEAR  CONSTRAINTS  ONLY.  NON¬ 
LINEAR  CONSTRAINTS  MAY  BE  USED  BY  CHANGING  THE  EQUATION  FOR 
V(S),  BUT  THEY  MUST  BE  CONVEX  IN  ORDER  FOR  THE  ALGORITHM  TO 
WORK . 
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(' 

C 

C 

Ctti^***t**tt*******1(**tt.'»t*M***t*****t*1^t***t***t****t***********r.***** 

VARIAHLE  DEFINITIONS 


REAL  VAR I ABIES: 

XCOORD  ^  I  hi:  X  COORraNATE  OF  THF  POINT  BEING  TESTED 
YCOORVi---  THE  Y -L'OURDINATE  OF  THE  POINT  BEING  TESTED 
ZCOORM-  THE  Z- COORDINATE  OF  THE  POINT  BEING  TESTED 
A(I)-  THE  COEFFICIENT  OF  X  IN  CONSTRAINT  I 
D(I)=  THE  COEFFICTENT  OF  Y  IN  CONSTRAINT  I 
C(I>^  THE  COEFFICIFNT  OF  Z  IN  CONSTRAINT  I 
K<I)^  THE  VALUE  OF  THE  RIGHT  HAND  SIDE  IN  CONSTRAINT  I 
V(I)=  THE  VALUE  OF  CONSTRAINT  I  EVALUATED  AT  THE  POINT 
BEING  TESTED 

INTEGER  variables: 

F=  THE  NUMBER  OF  FACES  (CONSTRAINTS)  ON  THE  REGION 
CHARACTER  VARIABLES: 

POINT=  THE  LOCATION  OF  THE  POINT?  EITHER  IN  OR  OUT 


SUPROUT I NE  CHEKPT  <  XCOORD  f  YCOORD  f ZCOORD , POI NT ) 

REAL  XCOORD, YCOORD, ZCOORD, A < 10), B( 10) ,C( 10) ,K<10> ,V( 10) 
INTEGER  F 
CHARACTER  P0INT*3 


—  IMPORTANT  NOTE  TO  USER  — 

THE  USER  MUST  SUPPLY  THE  NUMBER  OF  FACES  FOR  EACH  NEW  PROBLEM  HERE, 

Frr? 

THIS  LOOP  INITIALIZES  THF  COEFFICIENTS  AT  ZERO 
DO  500  S=1,F 
A(S)=0.0 
B<S)=0.0 
C ( S ) =0 . 0 
K(S)=0.0 
V<S)=0,0 
00  CONTINUE 

—  IMPORTANT  NOTE  TO  USER  — 

THF  USER  MUST  SUPPLY  THE  EQUATIONS  FOR  THE  CONSTRAINTS  IN  THE 
FORM  AX  +  BY  F  CZ  £  K  I.E,  HE  MUST  DETERMINE  AND  PROVIDE 
THE  VALUES  OF  A<S),  B(S),  CIS),  AND  K(S)  FOR  S*  1  TO  F,  FOR 
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r  jNST(^Nc:r:,  ir  iHf.  tiki'T  c:on'5traint  is  2X  i  .jy  +  4z  ^  5,  tul’n 
c  rtd)-:’,  K(l)-iy  f\ND  K(l)=5.  IF-'  THE  SECONIi  CONSTRAINT 

C  IS  X  1,  THEN  A(2)-l,  6(2)  ^0,  C(2)=0,  ANI.i  K(2)---l.  IF  A 

C  CONS  IRA]  NT  HAS  A  t3r<l:ATE.ft  THAN  SIGN,  THEN  REVERSE  THE  SIGN  BY 

C  NEGATING  THE  COEFFICIENTS.  FOR  INSTANCE,  X  +  Y  +  Z  >  3  BECOMES 

C  -X  -Y  -/  i,  SO  A<3)=-1,  B<3)=-1,  C<3)='l,  AND  K(3)=-l.  IF 

C  THE  CONSTRAINT  IS  AN  EQUALITY  CONSTRAINT,  THEN  FORM  TWO  EQUATIONS 

C  WITH  INEQUALITIES.  FOR  INSTANCE,  Y  T  Z  S  BECOMES  Y  +  Z  <  S 

r:  AND  -Y  -/  >  -s,  SO  A(4)=0,  B(4)=l,  C<4>  =  1,  K(4)=5,  A(fi)=0, 

C  D(5)  -I,  ('(b^^-1,  AND  K(S)=--5.  IF  YOU  HAVE  ANY  PROBLEMS, 

C  CONSULT  YOUR  LOCAL  MATH  EXPERT. 

C 

A  <  1  )  -  1 . 0 
B  <  2  >  ^  1 . 0 
n ( 3) ~ 1,0 
A<4)=-.l  .0 
B  <  5  )  =  - 1  . 0 
C<6)=-,1  .0 
A  <  7 )  ^  - 1 . 0 
B<7)=-.t  .0 
C(7)=-1.0 
K( 1 )=50.0 
K(2)=160.0 
K<3)*20.0 
K ( 4 ) »0 . 0 
K<5)=0.0 
K<6)=0.0 
K<7)=-2.3 

C  THE  CONSTRAINTS  AS  THEY  ARE  SHOWN  ABOVE  CORRESPOND  TO  CARTER  AND 
C  CARTER'S  CONCEPT  OF  THE  FEASIBLE  REGION  FOR  THE  COLOR  SPACING 
C  PROBLEM  AND  LOOK  LIKE  THIS  IN  RED,  GREEN,  AND  BLUE  COORDINATES: 

C  RED  ^  SO 

C  GF5EEN  <  160 

C  BLUE  20 

C  RED  >  0 

C  GREEN  2  0 

C  BLUE  2  0 

C  RED  +  GREEN  f  BLUE  2  2.3 

POINT='IN' 

C  THIS  LOOP  CHECKS  THE  POINT  AGAINST  EACH  OF  THE  CONSTRAINTS 
DO  510  G==I,F 

V  <  S ) = A  <  S ) *XCOORn+B  <  S ) #YCOORn+C  <  S ) *ZCOORD 
IF  <V(S)  .GT,  K(S))  THEN 
POINT=^'OUT' 

GO  TO  520 
END  IF 

510  CONTINUE 
520  END 
C 
C 


c  siiitfsuur INK  I  IX  PuiMi 

c 

c  (HF  purpose:  OE  this  SUltROUTTNF.  IS  TO  FIX  POINTS  ANH 

C  DETERMINE  WHETHER  OR  NOT  THE  PROGRAM  IS  DONE.  THE  LOGIC 

C  ASSOCIAIED  WITH  PiyiNG  POINT;'  IS  THAT,  ONCE  THE  MINIMUM 

C  DISTANCE  IS  MAXIMIZED,  THAT  DISTANCE  IS  NO  LONGER  UNDER  CONSIDER- 

C  ATION,  SO  ONE  OS  THE  ENDPOINTS  IS  FIXED  AI  ITS  LOCATION  AND  THE 

C  SEi.'OND  CLOSES!  PAIR  OF  POINTS  i;i  MAXIMIZED.  THIS  CONTINUES 

C  UNIIL  AIL  POINIS  ARE  FIXED,  AT  WHICH  TIME  THEY  ARE  UNFIXCD 

C  TO  SEE  ir  ANY  1 MPROOEMENTS  CAN  BE  MADE.  I! IE  PROGRAM  STOPS 

C  WHEN  NONE  CAN  BE  MADE. 

C 

C 

C*****i.M**t*****.***t*^^.**ttt****^**M******tti**t***t.t****************t* 

C 

C  UARIAPLE  DEFINITIONS 

C 

c  REAL,  oariadi.es: 

C  X(I)=  THE  X-COORDINATE  OF  POINT  I 

C  Y<I)=  the  Y-COORDINATE  OF  POINT  I 

C  Z<I)=  THE  Z-COORDINATE  OF  POINT  I 

C  DUMl  ,DUM2,DUM.J,DUM4=  DUMMY  OARIABLES  REQUIRED  FOR  SUBROUTINE  'FIND- 

C  DISTANCE'  BUT  NOT  USED  IN  THIS  SUBROUTINE 

C  MEHX<I)=  THE  X-COORDINATE  OF  POINT  I  FROM  THE  PREVIOUS  ITERATION 

C  MEMY<I)=  THE  Y-COORDINATE  OF  POINT  I  FROM  THE  PREVIOUS  ITERATION 

C  MEMZ(I)=  THE  Z-COORDINATE  OF  POINT  I  FROM  THE  PREVIOUS  ITERATION 

C  SIZE-  THE  STEP-SIZE 

C 

C  INTEGER  variables: 

C  N=  THE  NUMBFR  of  POINTS 

C  ENDPTl*  THE  NUMBER  ASSOCIATED  WITH  ENDPOINT  1 

C  ENDPT2--  THE  NUMBER  ASSOCIATED  WITH  ENDPOINT  2 

C  NUMFIX-  THE  NUMBER  OF  POINTS  THAT  HAVE  BFEN  FIXED 

C 

C  CHARACTER  VARIABLES*. 

C  STATPT<I)=  THE  STATUS  OF  POINT  IJ  EITHER  FIXED  OR  UNFIXED 

C  STATUS^  THE  STATUS  OF  THE  PROGRAM;  EITHER  DONE  OR  RUN 

C 
C 

c 

c 

SUBROUTINE  FIXPT ( X , Y , Z ,N , 5TATPT , ENDPTl ,ENDPT2 , 31 ZE , STATUS , MEMX , 
i  MEMYrMEMZ) 

REAL  X<50)  ,Y(50)  ,Z<50)  ,DUM1  ,DUM2<2)  ,  DUK3  <  2 )  ,  DLIM«»  (  2 )  ,  MEMX  ( 50 )  , 

S  MEMY<50) ,MEMZ<50) ,STZE 

INTEGER  N, ENDPTl, ENDPT2,NUMriX 
CHARACTER  STATPT ( 50) *1 ,STA TUS»4 
C 

C  —  IMPORTANT  NOTE  TO  USER  -- 

C 

C  IF  THE  USER  WISHES  TO  CHANGE  THE  MAXIMUM  STEP-SIZE,  HE  MUST  MAKE 
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o  n  n  o  n  n  n 


r  (iir.  L'Hi'SNnr  n'>  wi  i.  i  in  i:ir  i  kooham. 

c 

SX2E=110.0 

NLIMKIX^O 

C  THIS  LOOP  COUNIS  THr  NUMBER  OF  POINTS  THAT  ARE  ALREADY  FIXED 
DO  AOO  I---l,N 

ir  (STATPT<1)  .EQ.  'EM  THEN 
NnML'/X--NUMFTXH 
END  IF 

600  CUMITMUE 

CALL  I  IMDD  <  X  ,  Y  ,  Z  ,  M  ,  STATPT  ,  DUMl  ,  DUMZ  ,  DUM3 ,  DUM4  ,  E  NDPT  t  ,  ENDPTl’ , 

X  STAIUS) 

C  THESE  NEXT  TWO  IF  STATEMENTS  FIX  A  NEW  POINT 
IF  (STAtPT(ENDPTI)  .£Q.  'U')  THEN 
NUMFIX-NUMF IXTl 
STATPT  (C-:NDPT.I)-'F' 

GO  TO  610 
END  IF 

IF  <  STATPT  (E:NDPT2)  .EQ.  'U')  THEN 
NUMFIX=NUMFIXH 
STATPT<ENDPT2)-'F' 

END  IF 

C  IF  ALL  THE  POINTS  ARF  FIXED,  THEN  CHECK  TO  SEE  IF  THEY  ARE  IN  THE 
C  SAME  LOCATION  AS  BEFORE 
610  IF  (NUMFIX  .EQ.  N)  THEN 
DO  620  1=1, N 

IF  (MEMXd)  .NE.  X<I)  .OR.  MEMYCI)  .NE.  Y(I)  .OR. 

X  MEMZd)  .NE.  Zd>>  THEN 

STATUS*'RUN' 

GO  TO  630 

ELSE 

STATUS='DONE ' 

END  IF 

620  CONTINUE 

C  IF  THE  POINTS  HAVE  MOVED  SINCE  THE  LAST  ITERATION,  THEN  UNFIX  Al  I 
C  THE  POINTS,  STORE  THEIR  LOCATIONS  IN  MEMORY,  AND  RUN  THROUGH  THE 
C  PROGRAM  AGAIN 
630  no  640  1=1, N 

STATPTd  )=^'U' 

MEMXd)=Xd) 

MEMYd)aYd) 

MEMZd)=Zd) 

640  CONTINUE 

END  IF 
END 


********t*******t******$t***tM*tt******$***t**tt*t*****.*.*****t****M*M 

SUBROUTINE  PRINT 

THIS  SUBROUTINE  PRINTS  THE  LOCATION  OF  ALL  N  POINTS  TO  THE 
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('  K  [I.!.  I  .M  I  Lo  ;;i  i:i)ui '  i  tm  to  fAit,  rr  '  -  out  -  > . 

t: 

c 

C 

C  yORJABI.E  DEFINITIONS 

C 

C  IsFAl  VAI:I0I()  i-.s: 

r  xd'^  tiik:  x-roaRDiNAiF  of  i^oint  t 

C  Y<I>^  IMF  Y- COORD  (NrSTI:,  OF  POINI  1 

C  Z<I)  ^  THE  V  -COORDINrtTE  OF  POINT  I 

C 

c;  INTEGER  0  Al- :  I A  D L  E  S  : 

C  N-  THE  NUhDER  OF  POINTS 

C 

C 

C 
C 


710 

700 
720 
730 

C 
C 

c 

C  SUBROUTINE  RGBOUT 

C 

C  THIS  SUBROUTINE  CONVERTS  POINTS  FROM  THEIR  1.*,  U»,  AND  V* 

C  COORDINATES  TO  THEIR  TRISTIMULUS  VALUES  AND  THEN  TO  THEIR 

C  PHOSPHOR  LUMINANCES  AND  PRINTS  THESE  VALUES  TO  THE  FILE  NAMED 

C  'SPCOUT'  <A.K.A.  SPACE-OUT).  THE  USER  MAY  CHOOSE  TO  ESTABLISH 

C  DIFFERENT  VALUES  OF  MAXIMUM  LUMINANCE,  U-PRIME  NOT,  V-PRIME  NOT, 

C  AND  CHROMATICITY  COORDINATES  FOR  THE  GUNS  DEPENDING  ON  HIS  OR 

C  HER  PARTICULAR  APPLICATION.  ALSO,  FOR  SPACING  IN  A  SYSTEM 

C  OTHER  THAN  THE  CIE  L»U»V»  SYSTEM,  DIFFERENT  TRANSFORMATION 

C  EQUATIONS  MUST  BE  USED. 

C 

C 

Ctt****f.***t1lf1i*tt**tM******1li******M*t************t*****t1f**1f****tt1it** 

C 

C  VARIABLE  DEFINITIONS 

C 


SUBROUTINE  PRINT < X , Y , Z , N ) 

REAL  X(50) ,Y(50) ,Z<50) 

INTEGER  N 
WRITE(3,700) 

DO  710  T=t,N 

WrsITE<3,720)I,X(I)  ,Y<I)  ,Z(I) 

CONTINUE 

WRITE<3,730) 

FORMAT ( '  POINT' ,9X, 'L*' ,13X,'U*',13X,'V*'  ) 
r-0RMAT<3X,I2,7X,F8.3,7X,F8.3,7X,F8.3) 

FORMAT ('  ') 

END 
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ooonnnnnnnnonoonoonoonnnononnnn' 


c 

l-ilt''  VAI¬ 

Al'.l  1  1-  '  .  ■ 

r 

LS  rAp(  I  ' 

1  HP  vALur:  UP 

1  *  1  UK  i  OIMI  1 

c 

USTAI':<  I  ) 

rUL  VALUL  OP 

U*  PMR  POINT  I 

C 

VS  TAPs  (  1  > 

= 

THP  VALUE  OF 

V*  FOR  POINT  1 

c 

XR-  rHL 

X 

CHPOMAT IP  1 TY 

COORD I NA IE 

FOR 

1  PIP 

Rl  D  EiUN 

c 

YR--  IME 

Y 

CHPOMATIC  C lY 

COORD INA  IE 

FOR 

THE 

R:  D  GUN 

c 

XG=  I  UP 

X 

CinUlMAIlCI  TY 

COORD 1 NA IP 

FOR 

IHL 

GREEN  GUN 

c 

Yfi-  THL 

t' 

CHRnMATICTTY 

COORDTNA IE 

FOR 

I  HF 

ORELN  GUN 

c 

XL-  HIP 

X 

CHROMA  I ICI l  Y 

COORD (NA IP 

FOR 

I  Hi: 

BLUE  (lON 

c 

YJ<-  IHt 

f' 

CHROMA  I  It.'j:  )  Y 

COOh-MTNAIF 

;  OR 

IHE 

BLUE  GUN 

c 

YNGT^  THfe 

MAX  IMUM  1.  UMLNAMCL  ALI  OWAHLE 

yPHNOl  ■  THH-  UAl  LIE  Of  V -CRIME  MOT 
t.JCMNUT-=  THE  OMI  HE  01-  U  CK.1MI:  NOI 

XTOTftL([)^  THE  X--EOORliIMrtTE  OF  THE  TRISTIMULUG  OALIJE  FOR  POINT  I 

YTOTAL(J)-  iHE  Y  -  COORn  CNATE  OP  THP:  TRIGTIMULUS  VALUE  FOR  POINT  I 

7T0TAL<t)=  THE  Z  -COOKn INATE  OF  THE  TRISTIMULUS  VALUE  FOR  POINT  I 

l.)PftIMF(  I  )-  THE  VALUE  OP  IJ-PKIME  FOR  POINT  I 

VPRIMEd)^  THE  VALUE  OF  V--PRIME  FOR  POINT  I 

SMALLX(1)=  THE  X  CHROMATICITY  COORDINATE  OF  POINT  I 

SMALLY(I)=  THE  Y  CHROMATICITY  COORDINATE  OF  POINT  I 

SMALLZ<I)=  THE  Z  CHROMATICITY  COORDINATE  Of  POINT  I 

YRED(I)-  THE  RED  PHOSPHOR  LUMINANCE  FOR  POINT  I 

YGREEN<I)=  THE  GREEN  PHOSPHOR  LUMINANCE  FOR  POINT  I 

YBLUE<I)=  THE  BLUE  PHOSPHOR  LUMINANCE  FOR  POINT  I 

Kl=  THE  FIRST-ROW,  FIRST-COLUMN  VALUE  OF  THE  TRANSFORMATION  MATRIX 
THAT  CONVERTS  TRISTIMULUS  VALUES  TO  PHOSPHOR  LUMINANCES 
K2=  THE  FIRST-ROW,  SECOND -COLUMN  VALUE  OF  THE  MATRIX 
K3=  THE  FIRST-ROW,  THIRD-COLUMN  VALUE  OF  THE  MATRIX 
K4==  THE  THIRD-ROW,  FIRST-COLUMN  VALUE  OF  THE  MATRIX 
KS=  THE  THIRD-ROW,  SECOND-COLUMN  VALUE  OF  THE  MATRIX 
K6=  THE  THIRD-ROW,  THIRD-COLUMN  VALUE  OF  THE  MATRIX 
DET=  THE  DETERMINANT  OF  THE  TRANSFORMATION  MATRIX 

INTEGER  variables; 

N=  THE  NUMBER  OF  POINTS 

CHARACTER  VARIABLES: 

STATUS^  THE  STATUS  OF  THE  PROGRAM;  EITHER  'DUNE'  OR  'RUN' 


SUBROU  r I NE  RGBOUT ( LSTAR ,USTAR , VST AR , N , YRED , YGRFEN , YBLUE , STATUS ) 
REAL  LSTARISO) ,USTAR(50) ,VSTAR(50) ,XR,YR,XG,YG,XB,YB,YNOT, 

%  VPMNOT,UPMNOT,X  rOTAL  ( 50 )  ,  YTOTAL  (  50 )  , ZTOt^AL  ( 50 ) , UPRIME  ( 50 >  , 

S  VPRIME  <  50 ) , SMALLX  <  50 ) , SMALL Y ( 50 ) , SMALLZ ( 50 ) , YRED ( 50 ) , 

%  YGREEN(50) ,YBLUE<50) ,K1 , K2 , K3 , KA ,K5 ,K6 , DE T 

INTEGER  N 
CHARACTER  STATUS»4 
XR=.6 
YR-.3 


f,A 


.  .3 

YG-=.5B 

XB=  .  13 

Yh--  .OP 

YNOr -230.0 

IJPMNOT  :■  .  17f.37A 

OPMNOT  -  .  '^69/00 

K1 ^XR/YR 

K2-Xn/YG 

K3^^XP/YEi 

K4^( 1  -XR-YR)/YR 

K5:^<l-yG  YG)/YG 

K6=( 1-XR-YB)/YB 

riE  i=Kl»K6-Kl)KK5+K2*K4  -K2»K6+K3«K5-K3*K4 

DO  800  I-1,N 

C  THESE  ARE  THF  TRANSFORMATION  F.QUATIONS 

YTOTAL< I )=< ( <LSTAR< I )+16)/116)**3)»YN0T 

UPRIME<I)=(USTAR(I)/<13*LSTAR(I) ) ) hUPMNOT 

VPRIME ( I )  =  ( VSTAR < 1 ) / ( 13*LSTAR ( I ) >  >  +0PMN07 

SMALLX  < I ) =9*UPRIME  < I) / ( &*UPRIME ( I ) -1 6*VPRIME ( I )  + 1 2 ) 

SMALL Y  < I ) =4*VPR I ME ( I ) / ( 6»UPRIME  < I) - 1 6*VPR I ME (I )  + 1 2  > 

SMALLZ  < I )  =  1 -SMALLY  < I ) -SMALLX  < I > 

XTOTAL  < I ) aSMALLX ( I ) *YTOTAL  < 1 ) /SMALLY ( I ) 

ZTOTAL  < I ) =SMALLZ  < I ) *YTOTAL  < I > /SMALLY ( I ) 
YRFD<I)=<XT0TAL<I)*(K6-K5)+YT0TAL(I)*(K3»K5-K2»K&)+ 

S  ZT0TAL<I)*<K2-K3))/DET 

YGREEN<I)=<XT0TAL<I)*<K4-K&)+YT0TAL(I)*(K1»K&-K3»K4)+ 

I  ZT0TAL<I)*(K3-K1) )/DET 

YBLUE<I)*<XT0TAL<I)*<K5-K4)+YT0TAL(I)»(K2#K4-K1»K5)+ 

S  ZT0TAL<I)*<K1-K2))/DET 

800  CONTINUE 

C  THIS  SECTION  FORMATS  AND  PRINTS  THE  COMPUTED  TRISTIMULUS  VALUES  AND 
C  PHOSPHOR  LUMINANCES 

IF  (STATUS  .EQ.  'DONE')  THEN 
NRITE(3,810) 

WRITE<3,820) 

DO  830  1=1, N 

WRITE<3,840)I,XT0TAL<I),YT0TAL(I) rZTOTAL(I) 

830  CONTINUE 

WRITE<3,850) 

WRITE(3,860) 

WRITE<3,870) 

DO  880  1=1, N 

WRITE(3,840)I,YRED(I) ,YGREEN(I) ,YBLUE(I) 

880  CONTINUE 

WRITE<3,850) 

URITE<3,850) 

END  IF 

810  F0RMAT(14X,'THE  TRISTIMULUS  VALUES  ARE') 

820  F0RMAT(2X, 'POINT' ,9X, 'X' ,14X; 'Y' ,14X,'Z' > 

840  F0RMAT<3X,I2,7X,F8.3,7X,F8.3,7X,F8.3) 

850  FORMAT ('  ') 


^5 


860  r  Or'iNfU  '.?/,  ‘THfe.  I'HOSI-HOK-  (  UMlNrtNrfH  FDI:  THt  GUNS  rtHL 
870  ri)RKAT(7X, 'POINT' ,8X,  RED'  ,  1 1 X  ,  '  GtsEFN '  ,11X,'&L1)E'  ) 
END 
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Appendix  C 

EXAMPLE  OF  ROLEY  HEURISTIC 
SPACING  FOUR  POINTS  IN  A  SQUARE 


POINT  X-COORDINATE  Y-COORDINATE  Z-COORDINATE 


1 

.350 

.304 

.000 

2 

.  225 

.438 

.000 

3 

.118 

.247 

.000 

4 

.766 

.353 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.183 

CORRESPONDING  TO 

POINTS 

1  AND  2 

THE 

STEP-SIZE  IS 

.5000 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.850 

.804 

.  000 

2 

.225 

.938 

.000 

3 

.118 

.247 

.000 

4 

.766 

.353 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.458 

CORRESPONDING  TO 

POINTS 

1  AND  4 

THE 

STEP-SIZE  IS 

.5000 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.850 

.304 

.000 

2 

.225 

.'938 

.000 

3 

.118 

.247 

.000 

4 

.766 

.853 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.548 

CORRESPONDING  TO 

POINTS 

2  AND  4 

THE 

STEP-SIZE  IS 

.5000 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.850 

.304 

.000 

2 

.  225 

.938 

.000 

3 

.118 

.247 

.000 

4 

.766 

.853 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.548 

CORRESPONDING  TO 

POINTS 

2  AND  4 

THE 

STEP-SIZE  IS 

.2500 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.850 

.304 

.000 

2 

.225 

.938 

.000 

3 

.118 

.247 

.  000 

4 

.766 

.853 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.548 

CORRESPONDING  TO 

POINTS 

2  AND  4 

THE 

STEP-SIZE  IS 

.1250 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.  850 

.  304 

.  000 

2 

.100 

.938 

.000 
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3 

.118 

.247 

.000 

4 

.891 

.978 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.676 

CORRESPONDING  TO 

POINTS 

1  AND  4 

THE 

STEP-SIZE  IS 

.1250 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.  179 

.000 

2 

.100 

.938 

.000 

3 

.  118 

.247 

.000 

4 

.891 

.978 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.691 

CORRESPONDING  TO 

POINTS 

2  AND  3 

THE 

STEP-SIZE  IS 

.1250 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.  179 

.000 

2 

.  100 

.938 

.000 

3 

.118 

.122 

.000 

4 

.891 

.978 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.792 

CORRESPONDING  TO 

POINTS 

2  AND  4 

THE 

STEP-SIZE  IS 

.1250 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.179 

.000 

2 

.100 

.938 

.000 

3 

.118 

.122 

.000 

4 

.891 

.978 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.792 

CORRESPONDING  TO 

POINTS 

2  AND  4 

THE 

STEP-SIZE  IS 

.0625 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.179 

.000 

2 

.037 

.938 

.000 

3 

.118 

.122 

.000 

4 

.891 

.978 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.804 

CORRESPONDING  TO 

POINTS 

1  AND  4 

THE 

STEP-SIZE  IS 

.0625 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.116 

.000 

2 

.037 

.938 

.000 

3 

.118 

.122 

.000 

4 

.891 

.978 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.820 

CORRESPONDING  TO 

POINTS 

2  AND  3 

THE 

STEP-SIZE  IS 

.0625 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.116 

.000 

2 

.037 

.938 

.000 

3 

.118 

.059 

.000 

4 

.891 

.978 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.855 

CORRESPONDING  TO 

POINTS 

2  AND  4 

THE 

STEP-SIZE  IS 

.0625 
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POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.116 

.000 

2 

.037 

.938 

.000 

3 

.  118 

.059 

.000 

4 

.954 

.978 

.  000 

THE 

MAXIMIN  DISTANCE  IS 

.858 

CORRESPONDING  TO 

POINTS 

1  AND  3 

THE 

STEP-SIZE  IS 

.0625 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.054 

.000 

2 

.037 

.938 

.000 

3 

.056 

.059 

.  000 

4 

.954 

.978 

.  000 

THE 

MAXIMIN  DISTANCE  IS 

.878 

CORRESPONDING  TO 

POINTS 

2  AND  3 

THE 

STEP-SIZE  IS 

.0625 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.054 

.000 

2 

.037 

.938 

.000 

3 

.056 

.059 

.000 

4 

.954 

.978 

.  000 

THE 

MAXIMIN  DISTANCE  IS 

.878 

CORRESPONDING  TO 

POINTS 

2  AND  3 

THE 

STEP-SIZE  IS 

.0312 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.054 

.000 

2 

.006 

.969 

.000 

3 

.025 

.028 

.  000 

4 

.954 

.978 

.  000 

THE 

MAXIMIN  DISTANCE  IS 

.925 

CORRESPONDING  TO 

POINTS 

1  AND  4 

THE 

STEP-SIZE  IS 

.0312 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.023 

.  000 

2 

.006 

.969 

.000 

3 

.025 

.028 

.000 

4 

.954 

.978 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.941 

CORRESPONDING  TO 

POINTS 

2  AND  3 

THE 

STEP-SIZE  IS 

.0312 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.023 

.  000 

2 

.006 

.969 

.000 

3 

.025 

.028 

.000 

4 

.954 

.978 

•  000 

THE 

MAXIMIN  DISTANCE  IS 

.941 

CORRESPONDING  TO 

POINTS 

2  AND  3 

THE 

STEP-SIZE  IS 

.0156 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.023 

.  000 

2 

.006 

.969 

.000 

3 

.025 

.013 

.000 

4 

.954 

.978 

.000 
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THE 

MAXIMIN  DISTANCE  IS 

.948 

CORRESPONDING  TO 

POINTS 

2  AND  4 

THE 

STEP-SIZE  IS 

.0156 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.975 

.023 

.  000 

2 

.006 

.969 

.000 

3 

.025 

.013 

.000 

4 

.  969 

.994 

.  000 

THE 

MAXIMIN  DISTANCE  IS 

.950 

CORRESPONDING  TO 

POINTS 

1  AND  3 

THE 

STEP-SIZE  IS 

.0156 

POINT 

X-COORDINATE 

y-COORDINATE 

Z-COORDINATE 

1 

.991 

.023 

.000 

2 

.006 

.969 

.000 

3 

.025 

.013 

.000 

4 

.969 

.994 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.956 

CORRESPONDING  TO 

POINTS 

2  AND  3 

THE 

STEP-SIZE  IS 

.0156 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.991 

.023 

.  000 

2 

.006 

.984 

.000 

3 

.025 

.013 

.000 

4 

.969 

.994 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.963 

CORRESPONDING  TO 

POINTS 

2  AND  4 

THE 

STEP-SIZE  IS 

.0156 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.991 

.023 

.000 

2 

.006 

.984 

.000 

3 

.025 

.  013 

.000 

4 

.985 

.994 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.966 

CORRESPONDING  TO 

POINTS 

1  AND  3 

THE 

STEP-SIZE  IS 

.0156 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.991 

.007 

.  000 

2 

.006 

.984 

.000 

3 

.009 

.013 

.000 

4 

.985 

.994 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.972 

CORRESPONDING  TO 

POINTS 

2  AND  3 

THE 

STEP-SIZE  IS 

.0156 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATE 

1 

.991 

.007 

.000 

2 

.006 

.984 

.000 

3 

.009 

.013 

.  000 

4 

.985 

.994 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.972 

CORRESPONDING  TO 

POINTS 

2  AND  3 

THE 

STEP-SIZE  IS 

.0078 

POINT 

X-COORDINATE 

Y-COORDINATE 

Z-COORDINATF 

1 

.991 

.007 

.  000 
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2  .006 

.984 

.000 

3  .001 

.005 

.000 

4  .985 

.994 

.  000 

THE  MAXIMIN  DISTANCE  IS 

.979 

CORRESPONDING  TO  POINTS 

2  AND  4 

THE  STEP-SIZE  IS 
POINT  X-COORDINATE 


.0078 

y-COORDINATE 


Z-COORDINATE 


1  .991 

.007 

.000 

2  .006 

.992 

.000 

3  .001 

.005 

.000 

4  .993 

.994 

.000 

THE  MAXIMIN  DISTANCE  IS 

.987 

CORRESPONDING  TO  POINTS 

2  AND  4 

THE  STEP-SIZE  IS 


POINT 


X-COORDINATE 


.0078 

Y-COORDINATE 


Z-COORDINATE 


1  .991 

.007 

.  000 

2  .006 

.992 

.000 

3  .001 

.005 

.000 

4  .993 

.994 

.000 

THE  MAXIMIN  DISTANCE  IS 

.987 

CORRESPONDING  TO  POINTS 

2  AND  4 

THE  STEP-SIZE  IS 
POINT  X-COORDINATE 


.0039 

Y-COORDINATE 


Z-COORDINATE 


1  .991 

.007 

.000 

2  .002 

.996 

.000 

3  .001 

.005 

.000 

4  .997 

.998 

.000 

THE  MAXIMIN  DISTANCE  IS 

.989 

CORRESPONDING  TO  POINTS 

1  AND  3 

THE  STEP-SIZE  IS 
POINT  X-COORDINATE 


.0039 

Y-COORDINATE 


Z-COORDINATE 


1  .994 

.003 

.000 

2  .002 

.996 

.000 

3  .001 

.001 

.000 

4  .997 

.998 

.000 

THE  MAXIMIN  DISTANCE  IS 

.993 

CORRESPONDING  TO  POINTS 

1  AND  3 

THE  STEP-SIZE  IS 
POINT  X-COORDINATE 


.0039 

y-COORDINATE 


Z-COORDINATE 


1 

.998 

.003 

.000 

2 

.002 

.996 

.000 

3 

.001 

.001 

.  000 

4 

.997 

.998 

.000 

THE 

MAXIMIN  DISTANCE 

IS 

.994 

CORRESPONDING  TO  POINTS 

2  AND  4 

THE 

STEP-SIZE  IS 

.0039 

POINT 

X-COORDINATE 

y-COORDINATE 

Z-COORDINA' 

1 

.998 

.003 

.000 

2 

.002 

.996 

.000 

3 

.001 

.001 

.000 

4 

.997 

.998 

.000 

THE 

MAXIMIN  DISTANCE 

IS 

.994 

CORRESPONDING  TO  POINTS 

2  AND  4 
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THE 

STEP-SIZE  IS  .0020 

POINT 

X-COORDINATE  Y-COORDINATE 

Z-COORDINATE 

1 

.998 

.003 

.000 

2 

.000 

.998 

.000 

3 

.001 

.001 

.000 

4 

.999 

1.000 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.997 

CORRESPONDING  TO  POINTS 

1  AND  4 

THE 

STEP-SIZE  IS  .0020 

POINT 

X-COORDINATE  Y-COORDINATE 

Z-COORDINATE 

1 

.998 

.001 

.000 

2 

.000 

.998 

.000 

3 

.001 

.001 

.000 

4 

.999 

1.000 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.997 

CORRESPONDING  TO  POINTS 

2  AND  3 

THE 

STEP-SIZE  IS  .0020 

POINT 

X-COORDINATE  Y-COORDINATE 

Z-COORDINATE 

1 

.998 

.001 

.000 

2 

.000 

.998 

.000 

3 

.001 

.001 

.000 

4 

.999 

1.000 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.997 

CORRESPONDING  TO  POINTS 

2  AND  3 

THE 

STEP-SIZE  IS  .0010 

POINT 

X-COORDINATE  Y-COORDINATE 

Z-COORDINATE 

1 

.998 

.001 

.000 

2 

.000 

.999 

.000 

3 

.000 

.001 

.000 

4 

.999 

1.000 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.998 

CORRESPONDING  TO  POINTS 

2  AND  3 

THE 

STEP-SIZE  IS  .0010 

POINT 

X-COORDINATE  Y-COORDINATE 

Z-COORDINATE 

1 

.998 

.001 

.000 

2 

.000 

.999 

.000 

3 

.000 

.001 

.000 

4 

.999 

1.000 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.998 

CORRESPONDING  TO  POINTS 

2  AND  3 

THE 

STEP-SIZE  IS  .0005 

POINT 

X-COORDINATE  Y-COORDINATE 

Z-COORDINATE 

1 

.998 

.001 

.000 

2 

.  000 

1.000 

.000 

3 

.  000 

.000 

.000 

4 

.999 

1.000 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.998 

CORRESPONDING  TO  POINTS 

1  AND  3 

THE 

STEP-SIZE  IS  .0005 

POINT 

X-COORDINATE  Y-COORDINATE 

Z-COORDINATE 

1 

.999 

.001 

.000 

2 

.000 

1.000 

.000 

3 

.000 

.000 

.000 
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1.000 


000 


4  .999 

THE  MAXIMIN  DISTANCE  IS  .998 

CORRESPONDING  TO  POINTS  2  AND  4 
THE  STEP-SIZE  IS  .0005 

POINT  X-COORDINATE  Y-COORDINATE  Z-COORDINATE 


1  .999 

.001 

.000 

2  .000 

1.000 

.000 

3  .000 

.000 

.000 

4  .999 

1.000 

.000 

THE  MAXIMIN  DISTANCE  IS 

.999 

CORRESPONDING  TO  POINTS 

1  AND  3 

THE  STEP-SIZE  IS  .0005 

POINT  X-COORDINATE  Y-COORDINATE  Z-COORDINATE 


1  .999 

.001 

.000 

2  .000 

1.000 

.000 

3  .000 

.000 

.000 

4  .999 

1.000 

.000 

THE  MAXIMIN  DISTANCE  IS 

.999 

CORRESPONDING  TO  POINTS 

2  AND  4 

THE  STEP-SIZE  IS  .0005 

POINT  X-COORDINATE  Y-COORDINATE  Z-COORDINATE 


i  .999 

.001 

.000 

2  .000 

1.000 

.000 

3  .000 

.000 

.000 

4  1.000 

1.000 

.000 

THE  MAXIMIN  DISTANCE  IS 

.999 

CORRESPONDING  TO  POINTS 

1  AND  4 

THE  STEP-SIZE  IS  .0005 

POINT  X-COORDINATE  Y-COORDINATE  Z-COORDINATE 


1  1.000 

.001 

.000 

2  .000 

1.000 

.000 

3  .000 

.000 

.000 

4  1.000 

1.000 

.000 

THE  MAXIMIN  DISTANCE  IS 

.999 

CORRESPONDING  TO  POINTS 

2  AND  3 

THE  STEP-SIZE  IS  .0005 

POINT  X-COORDINATE  Y-COORDINATE  Z-COORDINATE 


1  1.000 

.001 

.000 

2  .000 

1.000 

.000 

3  .001 

.000 

.000 

4  1.000 

1.000 

.000 

THE  MAXIMIN  DISTANCE  IS 

.999 

CORRESPONDING  TO  POINTS 

2  AND  3 

THE  STEP-SIZE  IS  .0005 

POINT  X-COORDINATE  Y-COORDIHATE  Z-COORDINATE 


1  1.000 

.001 

.000 

2  .000 

1.000 

.000 

3  .001 

.  000 

.000 

4  1.000 

1.000 

.000 

THE  MAXIMIN  DISTANCE  IS 

.999 

CORRESPONDING  TO  POINTS 

2  AND  3 

THE  STEP-SIZE  IS  .0002 

POINT  X-COORDINATE  Y-COORDINATE  Z-COORDINATE 
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1 

1.000 

.001 

.  000 

2 

.000 

1.000 

.000 

3 

.000 

.000 

.  000 

4 

1.000 

1.000 

.000 

THE 

MAXIMIN  DISTANCE  IS 

.999 

CORRESPONDING  TO  POINTS 

1  AND  3 

THE 

STEP-SIZE  IS  .0002 

POINT 

X-COORDINATE  Y-COORDINATE 

Z-COORDINATE 

1 

1.000 

.000 

.  000 

2 

.000 

1.000 

.000 

3 

.000 

.000 

.000 

4 

1.000 

1.000 

.000 

THE 

MAXIMIN  DISTANCE  IS 

1.000 

CORRESPONDING  TO  POINTS  2  AND  4 
THE  STEP-SIZE  IS  .0002 


Although  the  solution  is  now  optimal,  the  program  continues  until  the  step- 
size  is  less  than  .0001  at  which  time  the  step-size  is  reset  back  to  .5  and 
point  number  three  is  fixed  in  its  place,  allowing  the  algorithm  to  concen¬ 
trate  on  trying  to  optimize  the  other  three  points.  The  step-size  then 
gradually  decreases  back  to  .0001  and  another  point  is  fixed,  and  so  on, 
until  all  four  points  are  fixed.  At  that  time,  the  locations  of  the  points 
are  stored  in  memory  and  all  the  points  are  unfixed  to  see  if  any  more 
improvements  can  be  made.  The  program  then  goes  through  the  whole  process 
again,  reducing  then  resetting  the  step-size  and  fixing  the  points  until  all 
four  points  have  been  fixed.  Because  no  improvements  could  be  made,  the 
program  is  done  at  that  time. 
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